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ABSTRACT 

1—5 In this paper, we derive scaling relations between photometric observ- 

QQ able quantities and disk galaxy rotation velocity V^ot, or Tully-Fisher re- 

lations (TFRs). Our methodology is dictated by our purpose of obtaining 
purely photometric, minimal-scattcr estimators of Vrot applicable to large 
galaxy samples from imaging surveys. To achieve this goal, wo have con- 
structed a sample of 189 disk galaxies at redshifts z < Q.l with long-slit 
Ha spectroscopy from Pizagno ct al. (2007) and new observations. By con- 
struction, this sample is a fair subsamplc of a large, well-defined parent disk 
sample of ~ 170 000 galaxies selected from the Sloan Digital Sky Survey 
Q Data Release 7 (SDSS DR7). The optimal photometric estimator of VJ-ot we 

^ find is stellar mass from Bell et al. (2003), based on the linear com- 

bination of a luminosity and a colour. Assuming a Kroupa IMF, we find: 

,Cd^ log[V8o/(kms"^)] = (2.142 ± 0.004) + (0.278 ± O.OlO)[log(Af,/Af0) - 10.10], 

where Vgo is the rotation velocity measured at the radius i?80 containing 80 per 
cent of the i-band galaxy light. This relation has an intrinsic Gaussian scatter 
^ a = 0.036 ± 0.005 dcx and a measured scatter (Tmoas = 0.056 dex in logVgo- 

For a fixed IMF, we find that the dynamical-to-stellar mass ratios within -Rsoj 
(Af(jyn/M*)(i?8o), decrease from approximately 10 to 3, as stellar mass in- 
creases from M* » 10^ to IO^^Mq. At a fixed stellar mass, (Mdyn/Mi.)(i?8o) 
increases with disk size, so that it correlates more tightly with stellar surface 

\^ density than with stellar mass or disk size alone. We interpret the observed 

variation in (Mdyn/Af«)(i?8o) with disk size as a reflection of the fact that disk 

^-H size dictates the radius at which M^y^/M^ is measured, and consequently, the 

fraction of the dark matter "seen" by the gas at that radius. For the lowest 
^ M* galaxies, we find a positive correlation between TFR residuals and disk 

sizes, indicating that the total density profile is dominated by dark matter 

/\ on these scales. For the highest galaxies, we find instead a weak negative 

H correlation, indicating a larger contribution of stars to the total density pro- 

file. This change in the sense of the correlation (from positive to negative) 
is consistent with the decreasing trend in (Afdyn/Af*)(i?8o) with stellar mass. 
In future work, we will use these results to study disk galaxy formation and 
evolution, and perform a fair, statistical analysis of the dynamics and masses 
of a photometrically-selected sample of disk galaxies. 

Key words: Key words: galaxies: spiral - galaxies: kinematics and dynam- 
ics. 
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In the standard cold dark matter (CDM) based model 
for disk galaxy formation, gas cools out of a hot gaseous 
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halo, maintaining its specific angular momentum, and 
forms a disk at the center of the potential well of a dark 



matter (DM) halo | |Fall fc Efstathiou|[T980[ ) . Although 
the basic picture has long been in place, a variety of 
physical processes underlying galaxy formation (such as 
star formation, feedback, angular momentum transfer, 
and mergers) is still poorly understood. 

It is thought that scaling relations between halo 
properties (e.g., halo mass and circular velocity) trans- 
late into scaling relations between observed galaxy prop- 
erties (e.g., stellar mass and disk rotation velocity). In 
detail, however, what one expects for the final relations 
(including the slope, zero-point, and scatter) depends 
on still-uncertain aspects of galaxy formation. 

Therefore, observed scaling relations for disk galax- 
ies are expected to provide strong constraints on galaxy 
formation scenarios. In particular, the well-established 
relation between disk rotation velocity and galaxy lu- 



minosity (TuUy & Fisher 19771 is arguably surpris- 



ingly tight, and hence, deserves special attention. In this 
work, we refer to scaling relations between disk rotation 
velocity and various galaxy properties (such as luminos- 
ity and stellar mass) generally as TuUy-Fisher relations 
(TFRs) . The purpose of studying these scaling relations 
has changed over the years, and the philosophy of sam- 
ple selection and calibration has varied accordingly. 

Originally, TFRs were used to determine galaxy 
distances and measure deviations from the cosmic flow 
(e.g.,|Tully fc Fisher|1977| [Strauss fc Willick|1995l[Gia 



vanel i et al.|1997||Courteau et al.|2000||Springob et al 



2009 1 . Accordingly, galaxy samples were pruned on the 
basis of morphological type to minimize the scatter in 
the relation, and ensure the smallest possible error in 
magnitude and distance. 

More recently, TFRs have been derived from galaxy 
samples including a broad range of morphological types 
(e.g., Courteau et al. 2007 Pizagno et al. 2007 here- 



after, P07), for the purpose of testing galaxy formation 
and evolution models (e.g., Dalcanton et al. 1997 Mo 



et al. 1998 Dutton et al. 2007 Gnedin et al. 20071 



These TFR studies attempt to constrain the theoreti- 
cal models, e.g., by quantifying, rather than minimizing, 
the scatter in the TFR. 

A third application of TFRs is connecting disk 
galaxies with their dark matter haloes through dynam- 
, [Dutton et al.|2010||Trujillo-Gomez et al.|2010 ) 



(e.g. 



Because long-slit spectroscopy is more expensive than 
imaging, one can exploit existing large imaging datasets 
by using TFRs to provide photometric estimates of disk 
rotation velocities Kot. Unlike previous TFR studies, 
we tailor our sample selection and calibration explicitly 
for this purpose. We consider several different photo- 
metric quantities O — including absolute magnitudes, 
synthetic magnitudes, stellar mass estimates, and bary- 
onic mass estimates — to determine the one that yields 
the tightest relation with Vrot, and hence, is the optimal 
estimator of Kot ■ 

In a future paper, we will combine results from 
this optimal TFR with weak lensing measurements to 



constrain the ratio of rotation velocities at the opti- 
cal and virial radii, Viot/Vvir- To avoid sample selec- 
tion issues, we want to use the same sample, or — in 
practice — similarly-defined samples, for the two analy- 
ses. The ratio Viot/K-ir provides a direct measurement 
of the slope of the total (dark -I- baryonic) mass pro- 
flle within the virial radius. Measurement of K-ot/Vvir 
can thus constrain halo structure and indicate whether 
baryons have significantly modified the halo potential 



(|Seljak|2002 


Dutton et al.[[2007l. 


Seljak I2OO2 


\ combined early TFR and weak lens 



ing measurements and inferred Vrot/Vvir = 1.8 with a 
2a lower limit of 1.4, for L* late-type galaxies. They 
found this result to be consistent with the prediction of 
a model with adiabatic contraction of the dark matter 
halo due to baryonic infall as described in |Blumenthal| 



et al.[ ( |1986| . More recently, [Dutton et al.[ ( [2010[ ) com- 
bined the TFR from P07 with halo mass measurements 
from different techniques (including weak lensing and 
satellite kinematics) and found Vrot/Vvir — 1, for disk 
galaxies with stellar masses aim 
to improve upon these measurements by carefully com- 
bining results from our TFR sample with weak lensing 
measurements on a similarly-selected lens galaxy sam- 
ple, leading to a result that is not as severely affected by 
systematic effects due to differences in sample selection, 
among others. 

We take advantage of the large, homogeneous 



dataset from the Sloan Digital Sky Survey (SDSS; York 
et al.|2000 l with well-defined photometry and available 



fibre spectroscopy. For the stacked weak lensing analy- 
sis, we will use a lens sample of ~ 10^ disk galaxies at 
redshifts z < 0.1; this is large enough to obtain decent 
S/N measurements. For the TFR analysis, we use a fair 
subsample of the lens sample, spanning the region in 
the parameter space of luminosity, galaxy colour, and 
size that the lens sample occupies. Our TFR sample 
consists of a set of galaxies with existing rotation curve 
measurements from P07, augmented by a comparable 
number of galaxies for which we have obtained new ro- 
tation curve measurements. By construction, our TFR 
sample matches the stellar mass function of the parent 
disk sample, and extends out to the same maximum 
redshift of 0.1 (unlike that used in P07). 

In addition to determining the slope, zero-point, 
and scatter in the TFR, we also study the residuals from 
the TFR, and their correlation with other galaxy prop- 



erties, such as disk size. Courteau & Rix ( 1999 1 argued 



that the lack of correlation between TFR residuals and 
disk size indicates that stars do not dominate the poten- 
tial well in the optical regions of disk galaxies (but also 
see [Dutton et al.|2007[ ). Both [Courteau et al.| ( [2007[ ) and 
P07 confirmed the lack of residual correlations in their 
galaxy samples. In this work, we go one step further and 
investigate residual correlations within different bins in 
stellar mass. 

We also constrain the ratio of dynamical to stel- 
lar masses within the optical radius of disk galax- 
ies, (Mdyn/M,)(_Ropt). We measure Mdyn(-Ropt) directly 
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from the rotation velocity measurements, and infer M, 
from photometric estimates based on galaxy luminosity 
and colour ( |Bell et al.|2"003l ). 

This paper is organized as follows: in !|2] we de- 
scribe previously existing data used in this work from 
SDSS (!j2]T]) and from P07 measurements {^^. In ijs] 
we describe our sample selection. In the next three sec- 
tions, we describe our derivation of the photometric and 
kinematic quantities used in the TFRs. In !j4] we de- 
scribe the bulge-disk decomposition fits from which we 
obtain disk parameters, such as disk size and axis ratio. 
In S|5] we define the photometric quantities O that we 
consider as estimators of disk rotation velocity. In ^ 
we describe the steps in the derivation of the kinematic 
quantities V^ot from the long-slit spectroscopy observa- 
tions. In |6.5[ we check for various systematic effects 
that may be affecting these measurements. 

In the next three sections, we analyse trends in 
these data. In !|7] we describe the modelling, fitting, 
and interpreting the TFRs. In ^ we present calibrated 
TFRs for the different photometric quantities consid- 
ered. In ^ we present alternative fits to the TFRs. 
In §10| we study correlations between residuals from 
the TFR and galaxy properties, such as disk axis ratio, 
galaxy colour, and disk size. In §11[ we calculate the 
dynamical-to-stellar mass ratios within the optical ra- 
dius of the galaxies, and study its correlations with stel- 
lar mass, disk size, and stellar surface density. Finally, 
we summarize and discuss our main results in §12| 

Throughout, we adopt a cosmology with ilm ~ 0.3, 
D.A = 0.7, and h = H/{100 kms-^Mpc"^) = 0.7, and 
express all lengths in physical (not comoving) coordi- 
nates. 



2 PREVIOUSLY-EXISTING DATA 

We use publicly-available data from the SDSS, and aug- 
ment long-slit spectroscopy data obtained by P07 with 



a set of new observations (described in S6.ll. In [2.1 



we describe SDSS imaging and spectroscopy, and define 
the SDSS photometric quantities used throughout this 
work. Then, in §2.2| we describe the sample and long-slit 
spectroscopy data of P07. 



2.1 SDSS data 



The seventh SDSS data release (DR7; 'Abazajia n|2009' I 
marks the completion of the SDSS-II phase of the sur- 
vey and includes 11,663 deg^ of imaging data, and over 
9380 deg'^ of spectroscopy. The imaging was carried 
out by drift-scanning the sky in photometric conditions 



(Hogg et al. 2001 Ivezic et al. 20041, in five bands 



(ugriz) (Fukugita et al. 1996 Smith et al. 20021 us- 



ing a specially-designed wide-field camera ( Gunn et al 



19981. All of these data were processed by automated 
pipelines that detect and measure photometric prop- 
erties of sources, and astrometrically calibrate the data 
( |Lupton et al.|2001|[Pie"r et al.|2003[ [Tucker et al.|2006 1. 



Objects are targeted for spectroscopy using the 
imaging data ( Blanton et al.|2003 1 . Main galaxy sample 
targets are selected as described by Strauss et al. ( 2002 I. 



The Main galaxy sample target selection includes a |Pet-| 
(19761 apparent magnitude cut of rp — 17.77 



mag, with slight variation in this cut across the sur- 
vey area. Targets are observed with a 640-fibre spectro- 
graph on the same telescope ( Gunn et al.||200^ l. These 
spectra are obtained through 3 -diameter fibres. In this 
work, SDSS fibre spectra are used only for sample se- 



lection (f3.1l; the extraction of kinematic information 



from these spectra will be explored in future work. 

In this work, we selected galaxy samples from the 
SDSS DR7 NYU- Value Added Galaxy Catalog (NYU- 



VAGC; Blanton et al. 20051. The selection criteria for 



inclusion in the VAGC is very similar to that for the 
Main galaxy sample. The imaging reductions used is 



a recent version, v5_4 (described in Abazajian 20041. 
The photometric calibration used is "uber-calibration" , 
which utilizes overlaps of SDSS runs to obtain a more 
consistent large-scale photometric calibration of the sur- 
vey ( Padmanabhan et al.|2008 l. 

We derive photometric quantities based on mea- 
surements from the SDSS imaging pipeline. These mea- 
surements include Petrosian and model apparent magni- 
tudes, Petrosian half-light radii -R50, and isophotal axis 
ratios giso- We note that both R50 and giso are not cor- 
rected for seeing, giao = &25/ffl25 or ISO_B/ISD_A are 
measured at an isophote of 25 mag arcsec"'^. They pro- 
vide a measure of the galaxy shape in its outer region^ 
and therefore a good initial estimate of disk axis ratios. 
We will only use isophotal axis ratios for sample se- 



lection (c.f. f3.2| and use seeing-corrected estimates of 



disk axis ratios and sizes from two-dimensional bulge- 
disk decompositions (described in for the TFR anal- 
ysis, in particular, for deriving inclination corrections to 
the rotation velocities. 

Absolute magnitudes are based on Petrosian ap- 
parent magnitudes, defined in [Strauss et al. (|2002[ ), fol- 



lowing the original proposal by Petrosian (19761. The 



essential feature of Petrosian magnitudes is that in the 
absence of seeing, they measure a constant fraction of a 
galaxy's light regardless of distance (or size). 

Galaxy colours are based on model apparent mag- 
nitudes, described in Stoughton et al. (2002') and'Abaza-J 



jian 



(20041. These magnitudes arc based on fitting the 



two-dimensional PSF-convolved galaxy image with ei- 
ther a pure deVaucouleurs or a pure exponential profile, 
depending on which model has a higher likelihood based 
on the r-band galaxy image. The fits to the other bands 
use the same model parameters — size, axis ratio, and 



^ Tlie semimajor axis of the 25 mag arcsec ^ isophote 025 is 
around 2 to 3 times the effective radius of the galaxy l |Ryden[ 
[2004) 1. 

In contrast, axis ratios from SDSS model fits are corrected 
for seeing, but are more severely affected by the light from 
the bulge than isophotal axis ratios. 
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position angle — as in the r-band galaxy image, to get 
stable colours. 

We correct both absolute magnitudes and colours 
for Galactic extinction using the dust maps of ISchlegel] 
et al.| ( |1998J , internal extinction (as described in ^5.1[ ), 
and fc-corrections to 2 = using the kcorrect prod- 



Blanton fc Roweis| ( |2007| )p] We 



uct version v4_l_4 of 

denote absolute magnitudes and colours that have not 
been corrected for internal extinction with the super- 
script "NC", e.g., M^°, Mf^, [g - r)^"^ . We reserve 
symbols without superscripts for quantities that have 
been corrected for internal extinction. 

In addition to these photometric quantities from 
the SDSS photometric pipeline, we also use'Sersic ( 1968) 
indices provided with the NYU-VAGC, determined 
from radial profile fits performed by |Blanton et al.| 
(20051). 



2.2 P07 data 

P07 selected 234 target galaxies from the SDSS DR2 
Main galaxy sample with redshifts z < 0.05, absolute 
magnitudes —22 < M^*^ < —18.5, and r-band isophotal 
axis ratios smaller than 0.6. In addition, they imposed 
an absolute magnitude-dependent upper limit on red- 
shift, so that the galaxies in their sample have apparent 
half-light radii larger than 2". In this work, we use a 
different set of selection criteria to define a disk galaxy 



sample (see i 3.1 1 



P07 observations were carried out at the Calar Alto 
Observatory using the TWIN spectrograph mounted on 
the 3.5 m telescope with 1200 lines mm~^ grating, 6200- 
73OOA spectral coverage, and 1.5" slit width, and at 
the MDM Observatory using the CCDS spectrograph 
mounted on the 2.4 m Hiltner telescope with 600 lines 
mm~^ grating, 6500-6994A spectral coverage, and 2" 
slit width. Integration times per galaxy vary from 1200 
s for bright galaxies to three exposures of 1200 s for 
faint galaxies. We refer the reader to the original pa- 
per for further details about the long-slit spectroscopy 
observations and data reduction. 

P07 obtained usable Hq rotation curves for 162 out 
of their 234 target galaxies (69 per cent); the other 72 
targets had insufficient Hq line emission for obtaining 
rotation curves. Out of these 162 galaxies, we found 
that four were flagged for redshift measurement warn- 
ings, and one is not included in the NYU-VAGC. Out 
of the remaining 157 galaxies (hereafter referred to as 
the "P07 galaxy sample"), 99 pass our sample selec- 



tion criteria (described in S 3.1 1 and are included in our 
analysis. For these galaxies, we re-processed the long- 
slit spectroscopy data obtained by P07 with the same 
pipeline used for data from our new observations. 



3 SAMPLE SELECTION 

As motivated in |T] we define two disk galaxy samples- 
(i) a parent disk sample of 169 563 galaxies, large enough 
to allow stacked weak lensing analysis, as well as de- 
tailed statistical analysis of galaxy properties, and (ii) 
a child disk sample of 189 galaxies, a subsample of the 
parent disk sample with measured rotation curves for 
kinematic analysis. 

Insofar as this child disk sample is a fair subsample 
of the parent disk sample, results derived from it can 
be applied to the larger sample. Thus, we have made an 
effort to construct a child disk sample that spans the 
parameter space of galaxy properties — colour, physical 
size, and galaxy type — covered by the parent disk sam- 
ple. Our child disk sample is composed of 99 galaxies 



from the P07 galaxy sample (see [ 2.2 1 and 90 galaxies 



for which we acquired new long-slit spectra, totalling 
189 galaxies with Ha rotation curves. 

We describe the selection criteria used to define the 
parent disk sample in §3.1[ We describe the construction 
of the child disk sample in §3.2[ Finally, we show that 
the distributions of galaxy properties are the same for 
the two samples in §3.3[ 



3.1 Parent disk sample 

We define the parent disk sample to be a large sample of 
nearby, star-forming, disk galaxies from the SDSS, ap- 
propriate for weak lensing and detailed statistical anal- 
ysis. Galaxy colour and Sersic index are often used in 
the literature to select disk galaxies, but applying such 
cuts needlessly excludes a significant population of disk 
galaxies ( Mailer et al.|2009 l. Thus, we do not apply any 
cuts on galaxy colour, and apply only a conservative cut 
on Sersic index. 

To select star-forming galaxies, we primarily rely 
on the Ha emission-line strength observed through the 
SDSS fibre as a marker of star formationj^Thus, we ex- 
plicitly select galaxies that are likely to have sufficient 
Ha emission for obtaining rotation curves. As a con- 
sequence, gas-poor SOs will be excluded automatically. 
This selection enabled us to efficiently use our observing 
time and obtain usable rotation curves for all of our tar- 
gets (compared to only 69 per cent of targets for P07). 

We select the parent disk sample from the SDSS 
DR7 NYU-VAGC described in ^ To clean the sam- 
ple, we exclude candidate objects that have spectro- 
scopic redshift measurement warnings (0.8 per cent) and 
those that are not spectroscopically confirmed as having 
class "GALAXY" (2.5 per cent). This selection yields a 
total of 686 656 galaxies. 

We then apply the following selection criteria: 



^ In other words, galaxy colours g — r = M^ad,g ~ ^mod.n 
where M^ad,g ^ud M^^^ ,, are the g and r-band model ab- 
solute magnitudes corrected for both Galactic and internal 
extinction and fc-corrected to 2 = 0, respectively. 



* The 3" SDSS fibres probe the integrated light from the 
central regions of nearby galaxies (at redshifts of 0.02, 0.07, 
and 0.10, the fibre covers the central 1.2, 4.0, and 5.6 kpc, 
respectively). 
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(i) 0.02 < z < 0.10 

(ii) -22.5 < M^'^ < -18.0 

(iii) /ob^Ha) > 2 X W~^^ erg s-^cm-^ 

(iv) 0.5 < < 5.9 AND < 1.7 - {Mr + 18.0) 

/ N , / [QUI] 5007 \ , ajil I 1 1 Q 

(V) log 1^ H/3 y log([NII]6583/Ha)-0.47 + ^'^^^ 

where z is the galaxy redshift, M^*^ is the r-band ab- 
solute magnitude corrected for Galactic extinction and 
fc-corrected to z = (but not corrected for internal ex- 
tinction), /oba(Ha) is the observed Ha emission-line flux 
through the SDSS spectroscopic fibre, and Ua is the ri- 



band Sersic index determined by Blanton et al. ( 2005 1 



The final cut removes galaxies with active nuclei (AGN) 
and is based on the ratio of high- to low-ionization 
emission-line fluxes, as measured through the SDSS fi- 
bre. We describe each of these cuts in detail below. 

(i) Redshift cut. We apply a lower limit on red- 
shift that ensures that the uncertainty in the absolute 
magnitude due to peculiar velocities (~ 300 kms~^) is 
no larger than 51og(l -I- 300/6000) = 0.11 mag. We also 
note that galaxies with z < 0.02 are not very useful for 
weak lensing because they would have negligible tangen- 
tial shear compared to lens galaxies at higher redshifts. 

The maximum redshift of 0.1 is a practical choice 
given the angular resolution required to obtain resolved 
rotation curves, as well as the need to obtain long- 
slit spectroscopic observations for the child disk sam- 
ple within a reasonable time. After applying this cut, a 
total of 327 027 candidates remain. 

(ii) Absolute magnitude cut. Following P07, the 
bright-end cut at M^^ — —22.5 mag excludes many 
non-starforming (i.e., elliptical) galaxies that may con- 
taminate the sample. This cut removes 10 852 candi- 
dates or 3.3 per cent of those from (i). 

(iii) Ha flux cut. We apply a lower limit on the 
Ha emission-line flux observed through the SDSS spec- 
troscopic fibre. We measure the Ha flux, /obs(HQ), by 
simultaneously fitting a set of three Gaussians plus a lin- 
ear continuum to the Ha and [NH]6548,6583 lines (with 
the ratio of [NH]6548 to [NII]6583 flux fixed to 1:3 and 
all line widths required to be equal). 

We choose the lower limit to be the value that 
roughly divides the bimodal distribution in Ha flux 
formed by the non-starforming and star-forming galaxy 
populations. This cut removes 37 per cent of candidates 
from (ii), leaving a total of 198 769 galaxies. 

We note that this cut introduces some redshift de- 
pendence. Since for nearby galaxies, the SDSS fibre is 
only sensitive to the central regions, some star-forming 
galaxies with extended Ha emission but suppressed cen- 
tral Ha emission would be rejected. 

(iv) Sersic index cut. We apply a conservative 
upper limit on the i-band Sersic index that depends 
on absolute magnitude. This cut closely approximates 



the 90th percentile Ua values at each MJ!**^ (increasing 
linearly from 2.2 at —18.5 mag to 5.7 at —22.0 mag). We 
also restrict to galaxies with 0.5 < rig < 5.9 to ensure 
that the fits are not at the limits of the the allowed 
parameter space, Us G [0.5, 6.0]. This cut removes 10 019 
candidates or 5 per cent of those from (iii). 

(v) AGN removal cut. It is important to exclude 
active galaxies (AGN) from our sample to ensure that 
most or all of the observed Ha emission comes from the 
disk, rather than the central engine, and can thus reli- 
ably trace the rotation of the disk. We use a standard 



diagnostic emission-line ratio diagram (Baldwin et al. 



1981 1 and apply the theoretical emission-line ratio cut 
of Kewley et al. |2001 ) to classify AGN. This cut is more 
conservative (removes fewer galaxies) than the empirical 
cut of Kauffmann et al. (20031. In line with being con- 



servative, we only exclude galaxies that have S/N > 1 
on all six emission lines. This final cut removes 10 per 
cent of candidates from (iv), and yields our parent 
disk sample with 169 563 galaxies. 

To measure the emission-line fluxes, we first sub- 
tract a best-fit stellar continuum from the SDSS spec- 
trum, following the procedure of Hao et al. (20051. 
We then perform a Gaussian fit to the continuum- 
subtracted spectrum similar to the procedure described 
above to measure the Ha flux. We repeat the process 
to measure the Rp and [OIII]4959,5007 lines (with the 
ratio of [OIII]4959 to [OIII]5007 flux fixed to 1:3). 

We have not applied an axis ratio cut to the parent 
disk sample, but we require r-band isophotal axis ratios 
Qiso < 0.6 on the child disk sample, following P07, to 
allow accurate inclination corrections to the galaxy ro- 
tation velocities. This cut is satisfied by 45 per cent of 
the parent disk sample, or 75 668 galaxies. 

Figure [l] shows the distributions of galaxy proper- 
ties for the parent disk sample galaxies (all galaxies and 
those with giso < 0.6 are shown in thick and thin black 
histograms, respectively). Also shown are the distribu- 
tions for galaxies that satisfy the redshift cut (i) but 
were excluded by the rest of the selection criteria (red 
histograms). Panels show the following galaxy proper- 
ties: absolute magnitude Mr, stellar mass M«, Petrosian 
apparent magnitude rp, redshift z, g — r colour, physical 
i-band Petrosian half-light radius i?5o,phys, i-band Sersic 
index ng, and r-band isophotal axis ratio qiso- Stellar 



masses M* correspond to the Kroupa (2002| IMF and 
were determined from i-band absolute magnitudes and 
g — r colours, using stellar mass-to-light ratio estimates 



from Bell et al. { 2003 1, as described in i 5.4.1 The other 
photometric quantities were defined in pT) 

The distributions in g — r and rig indicate that our 
selection criteria have successfully excluded the popula- 
tion of non-starforming (elliptical) galaxies, which con- 
stitute the red peak in g — r colour at ~ 0.75 mag and 



^ This is a reasonable fraction of AGN for a sa mple of 
emission-line galaxies. Out of ~ 10^ SDSS galaxies, |Brinch-| 
|mann et al. (2004) classified 68.4 per cent as emission-line 
galaxies, and 10.4 per cent of those as AGN. 
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Figure 1. Distribution of basic properties of galaxies in tiic parent disk sample (thick black histograms; 169 563 galaxies), those 
in the parent disk sample with isophotal axis ratios q^^^ < 0-6 (thin black histograms; 75 668 galaxies), and those galaxies that 
satisfy the redshift cut (i) but were excluded from the parent disk sample by the rest of the selection criteria (red histograms; 
157 464 galaxies). Upper panels show absolute magnitudes M^, stellar masses M*, Petrosian apparent magnitudes rp, and 
redshifts z. Lower panels show g — r colours, physical Petrosian half-light radii iJ5o,physi Sersic indices ris, and isophotal axis 
ratios qiso- 



the broad bump of high Sersic indices around ris ~ 3.5 
(see red histograms in Fig. [l] for the excluded galaxies). 
Recall that our selection relies primarily on the presence 
of the Hq line emission and does not impose any strict 
cuts in either colour or Sersic index. 



3.2 Child disk sample 

We aim to construct a child disk sample that is a fair 
subsample of the parent disk sample for use in our kine- 
matic analysis. By combining existing data from P07 
with newly-acquired observations, we have constructed 
a child disk sample of 189 disk galaxies with usable Hq 
rotation curves. 

Of the 157 galaxies in the P07 galaxy sample, 99 



pass our selection criteria (in S3.ll and are included in 
the child disk sample. Of the 58 rejected galaxies, 3 
failed due to the absolute magnitude cut, 17 due to the 
Hci flux cut, 11 due to the Sersic index cut, and 27 due to 
the AGN cut. The AGN fraction of 17 per cent is larger 
than the overall fraction of 10 per cent (c.f. criterion 



(v) in S 3.1 1. This may be attributed to the tendency for 
galaxies in the P07 sample to be of earlier type, and 
therefore m ore hkely to host AGN ( |Bundy et aI1|2008| 
but also see Xue et al.|2010 ). 

The other 90 galaxies in the child disk sample have 
been selected from the parent disk sample to fill in re- 
gions of the parameter space in luminosity, colour, and 
size that were not well covered by the existing dataset 
(see Fie. |2|. We require these galaxies to have r-band 
isophotal axis ratios giso < 0-6, following P07, to avoid 
large inclination corrections to the rotation velocities. 
We have targeted and acquired long-slit observations of 
these galaxies, as described in §6.1[ 

Table [T] lists the 189 galaxies in the child disk sam- 
ple, together with their basic properties and the instru- 
ment used in their observation. For the five galaxies 
with repeat observations, instruments from both obser- 
vations are listed. 

Figure [2] shows the distribution of child disk sample 
galaxies in the parameter space of absolute magnitude 
Mr, g — r colour, and physical half-light radius i?5o,phys 
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Table 1. Basic properties of the 189 galaxies in the child disk sample. 





Galaxy name 


Spec. ID 


2 


Mr 


g-r 


^50,phys 


riB 


Qiso 


Instr. 










(mag) 


(mag) 


(kpc) 










(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


SDSS 


J001006.61-002609.7 


0389-303-51794 


0.0321 


-21.13(0.10) 


0.32(0.03) 


2.42(0.08) 


1.22 


0.42 


TWIN 


SDSS 


J001708.75-005728.9 


0390-300-51900 


0.0189 


-19.62(0.13) 


0.57(0.04) 


3.13(0.17) 


1.45 


0.57 


TWIN 


SDSS 


J002844.82+160058.8 


0417-329-51821 


0.0947 


-22.37(0.23) 


0.33(0.06) 


6.08(0.06) 


1.71 


0.17 


DIS 


SDSS 


J003112.09-002426.4 


0391-063-51782 


0.0194 


-20.41 (0.12) 


0.76 (0.03) 


2.00(0.10) 


2.95 


0.53 


TWIN 


SDSS 


J004916.23+154821.0 


0419-602-51879 


0.0846 


-20.31 (0.07) 


0.54 (0.04) 


5.65(0.07) 


1.37 


0.53 


DIS 


SDSS 


J004935. 71+010655. 2 


0394-380-51913 


0.0176 


-21.21 (0.26) 


0.82(0.03) 


5.40(0.31) 


2.42 


0.17 


TWIN 


SDSS 


J011750.26+133026.3 


0423-044-51821 


0.0326 


-20.01 (0.07) 


0.36(0.03) 


3.83(0.12) 


1.44 


0.53 


DIS 


SDSS 


J012317.00-005421.6 


0399-254-51817 


0.0259 


-20.92(0.10) 


0.62 (0.03) 


2.33(0.09) 


2.02 


0.39 


TWIN 


SDSS 


J012340.12+004056.4 


0399-478-51817 


0.0334 


-21.23(0.14) 


0.60 (0.03) 


3.01 (0.09) 


2.35 


0.32 


TWIN 


SDSS 


J012438.08-000346.5 


0399-178-51817 


0.0277 


-20.92(0.14) 


0.64 (0.03) 


3.61 (0.13) 


1.95 


0.35 


TWIN 



Notes. — Col. (1): SDSS name of galaxy. Col. (2): Spectroscopic ID: plate, fiber ID, and MJD of SDSS fibre spectrum. Col. 
(3): Redshift. Col. (4): r-band Petrosian absolute magnitude, fc-corrected to 2 = 0, and corrected for Galactic and internal 
extinction. Col. (5): g — r model colour, fc-corrected to 2 = 0, and corrected for Galactic and internal extinction. Col. (6): 
physical i-band Petrosian half-light radius. Col. (7): i-band Sersic index. Col. (8): r-band isophotal axis ratio. Col. (9): 
Instrument(s) used for observation- TWIN/CCDS (P07), DIS (this work). Quantities in parentheses are 1-a uncertainties. 
The uncertainty in ^?50,phys includes a distance uncertainty assuming a peculiar velocity of 300 kms~^. 




Figure 2. Distribution in the parameter space of absolute magnitude Mr, g — r colour (left panel) and physical half-light radius 
R 

50,phys (right panel) of 189 galaxies in the child disk sample. Newly-observed galaxies are shown as open circles (90 galaxies), 
and previously-observed galaxies from P07 are shown as filled circles (99 galaxies). Galaxies in the parent disk sample (169 563 
galaxies) are shown as yellow dots, and their median relations and 1- and 2-a percentile bounds are shown in solid, dashed, and 
dotted lines, respectively. For comparison, galaxies from the P07 galaxy sample that failed our selection criteria (58 galaxies) 
are shown as red stars. 



(open and filled circles), compared with that of parent 
disk sample galaxies (yellow dots). Also shown are the 
median relations (solid lines) and the 1- and 2-a per- 
centile limits (dashed and dotted lines, respectively) for 
the parent disk sample. As desired, the child disk sam- 
ple spans the parameter space covered by the parent 
disk sample fairly well over most of the luminosity range 
covered. Moreover, 68 and 93 per cent of galaxies in the 
child disk sample lie within the 1- and 2-a percentile 
bounds of the parent disk sample, close to the expected 
values of 68 and 95 per cent, respectively. 

In Fig. [2j previously-observed and newly-observed 



galaxies are shown as filled and open circles, respec- 
tively. The former cover the bright and faint ends, while 
the latter fill in the region of intermediate absolute 
magnitudes —21 < Mr < —19.5 mag, by construction. 
Galaxies from the P07 galaxy sample that failed our se- 
lection criteria are shown as red stars. These excluded 
galaxies tend to be more luminous and redder than the 
mean population (i.e., earlier-type disks). 

Figure [3] compares the distribution of galaxy prop- 
erties for the child, P07, and parent disk galaxy samples 
(thick black, thin black, and red histograms, respec- 
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Figure 3. Distribution of basic properties for galaxies in the child disk sample (thick black histograms; 189 galaxies), P07 
galaxy sample (thin black histograms; 157 galaxies), and parent disk sample (red histograms; 169 563 galaxies). As in Fig. [T] 
upper panels show absolute magnitudes Mr, stellar masses M*, Petrosian apparent magnitudes rp, and redshifts z. Lower 
panels show g — r colours, physical Petrosian half-light radii -f?5o,physi Sersic indices rtg, and isophotal axis ratios qiso- 



Table 2. Mean galaxy properties for the parent disk, child 
disk, and P07 galaxy samples. 



brighter luminosities, higher stellar masses, and higher 
Sersic indices (but not redder colours). 



Sample 


N 




(Mr) 




{Rso) 








(mag) 


(mag) 


(kpc) 


Parent 


169 563 


0.0637 


-20.23 


0.53 


3.35 


Child 


189 


0.0407 


-20.50 


0.53 


3.49 


P07 


157 


0.0289 


-20.70 


0.61 


3.43 



tively). Table |2] compares mean values of some basic 
galaxy properties for the three samples. 

Figures [2] and [S] show that the P07 galaxy sam- 
ple is skewed toward brighter luminosities, higher stellar 
masses, redder g — r colours, and higher Sersic indices, 
compared to the parent and child disk samples. By con- 
struction, the child disk sample follows the distribution 
of the parent disk sample more closely. However, be- 
cause the child disk sample includes many galaxies from 
the P07 galaxy sample, it retains a slight bias towards 



3.3 Child vs. parent disk sample 

We check that the child disk sample is a fair subsam- 
ple of the parent disk sample once the small difference 
in their luminosity distributions are taken into account. 
To do this, we construct ten random subsamples of the 
parent disk sample with the same number of galaxies 
and the same Mr distribution (in 0.3 mag wide bins) 
as the child disk sample. We find that the distributions 
in g — r colour, Sersic index Ua, and physical half light 
radius 7?5o,phys for the child disk sample are consistent 
with the mean distributions for the ten random sub- 
samples of the parent disk sample within the sampling 
variance (as can be gleaned from Fig. [3|. 

We also perform a Kolmogorov-Smirnov (KS) test 
to check whether each pair of distributions is consis- 
tent with being drawn from the same underlying pop- 
ulation. We find mean D-statistic values equal to 0.12, 
0.11, and 0.14 for the distributions in g — r, ris, and 
-Rso.phys, respectively {D = if the two distributions are 
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identical). The large values of the corresponding prob- 
abilities P{D > observed) — 15, 26, and 9 per cent, 
respectively — indicate that we cannot reject the null 
hypothesis that the two distributions were drawn from 
the same population. 

In contrast, the P07 galaxy sample is highly incon- 
sistent with being drawn from the same distribution as 
the parent disk sample. Performing the same matching 
in absolute magnitude distributions and Kolmogorov- 
Smirnov tests, we find D-statistic values equal to 0.20, 
0.27, and 0.24, all with probabilities P{D > observed) 
much less than 1 per cent — 7 x 10~^, 9 x 10~^, and 
8 X 10"*. 



4 BULGE-DISK DECOMPOSITIONS 

First, we describe our motivation for performing bulge- 
disk decomposition (B-l-D) fits and how these results are 
used in this work (^4.1 1. Then, we describe the B-l-D 
fitting procedure {[4.2 \. Finally, we show example fits 
and discuss some results (S4.3I. 



4.1 Motivation and use 

We perform two-dimensional B-l-D fits of galaxies in the 
child disk galaxy sample to determine their disk sizes 
and axis ratios, and enable accurate rotation velocity 
estimates. The fits also provide disk-to-total light ra- 
tios, D/T, and surface brightness profiles, which we will 
present and analyze in future work. 

In this work, we use only three basic parameters 
from the B-l-D fits: 

• Qd'. disk axis ratio 

• -Rd: disk scale length 

• Rso'- radius containing 80 per cent of the total 
galaxy light 

We fit to the SDSS i-band galaxy images because this 
band is less affected by dust than the g and r bands 
(while the z band has noisier photometry). 

Disk axis ratios from the B-l-D fits are more ac- 
curate than the isophotal axis ratios giao from the SDSS 
photometric pipeline (used in our sample selection), be- 
cause the latter is affected by the presence of a bulge 
and variations in seeing conditions. We therefore use qd 
to determine internal extinction corrections to absolute 
magnitudes and colours (as described in j j5.1| |. More im- 
portantly, we use qd to determine inclination corrections 



to the observed rotation velocities (c.f. Eq. 22 1. 



4.2 Fitting procedure 

Numerous studies have performed two-dimensional 
bulge-disk decomposition of galaxy light profiles, in or- 



der to quantify their morphological properties ( Gadotti 
20091 [Benson et al][2007l [Allen et al.][2U06l [Peng et al 
2002'). However, a bulge-disk decomposition into an ex- 
ponential disk plus a Sersic bulge requires a fit with a 



large number (12) of free parameters, many of which 
are degenerate, in particular, the bulge Sersic index, 
the bulge effective radius, and the bulge central surface 
brightness. 

Fortunately, the surface density profiles of nearby 
disk galaxies are adequately described, on average, by a 



double-exponential profile ([MacArthur et al. 


2003 Gra- 


ham & Worley|2008 


McDonald et al. 


20111 


^ Since we 



are not interested in the detailed morphology of galaxy 
bulges in this work, we have chosen to fix the bulge 
Sersic index to 1, yielding more robust fits to the bulge 
and diskQ 

We fit the i-band galaxy images with a double- 
exponential distribution: 

I{R) = Id exp(-J?/iid) + h exp(-1.68i?/i?eff), (1) 

where Rd is the disk scale length, R^s is the effective 
radius of the bulge (equal to 1.68 times its scale length), 
and Id, Ih are the central surface brightness of the disk 
and bulge components, respectively. The radial coordi- 
nate 7?, is given by an ellipse centered at {xq, yo) 



R^{{a 



- xo)^ + {{y - yo)/q(b,d)f)^''^ 



(2) 



where q{b,d) denotes the axis ratios of the bulge and disk, 
respectively. As with the SDSS exponential models, the 
surface brightness of each component is truncated at 4 
times the effective radius, with a smooth fall-off to zero 
between 3 and 4 effective radii. 

We use the SDSS DR7 i-band atlas images from 
the most recent reductions (v5_6) ( |Aihara et al.|20Tl] |. 
These are sky-subtracted imagesj of typically less than 
200 X 200 pixels. The fits are done by minimizing the 
weighted difference between the image and the model 
convolved with the PSF at the galaxy position. The 
minimizer used is the Levenberg-Marquardt minimizer 



mpfit2dfun in IDL (Markwardt 2009 1. Each pixel is 
weighted by its inverse variance, which is computed us- 
ing the locally measured sky background counts from 
SDSS for each galaxy. 

The initial conditions for the fit are taken from 
the SDSS measurements of R_EXP, the exponential scale 
length, and AB_EXP, the exponential fit axis ratio, and 
the total fiux of the atlas image. The initial axis ratio 
and position angle of the bulge are chosen to be the 
same as those of the disk. The initial scale length of the 
bulge is set to 40% of the disk, and the initial D/T is set 



^ From their B+D fits of nearby disk galaxies, [MacArthur[ 
[et al.[ | [2003[ l found bulge Sersic indices ranging from 0.1 to 
2, with a mean very close to 1. 

The degeneracies are not completely removed, but fixing 
the bulge Sersic index decreases the covariance between bulge 
flux normalization and bulge scale length. 
® Improper sky subtraction will afli'ect the B-l-D fits. The 
general trend in SDSS DR7 (and DR8) is for the sky to be 
overestimated near bright, large galaxies, which leads to un- 
derestimation of the galaxies' sizes and luminosities. These 
sky subtraction issues are not expected to affect the typical 
galaxy in our sample. 
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to 0.6. The fits are robust against changes in the initial 
conditions, specifically changes in the initial bulge size 
and D/T. 

In order to compare the model to the atlas im- 
age during fitting, the model is convolved with the lo- 
cally measured point-spread-function from the SDSS 
pipeline. For each galaxy, the fitter returns the best- 
fit parameters, the covariance matrix of these parame- 
ters, and the reduced value for the fit. Since each fit 
consists of two exponential profiles, we chose the expo- 
nential with the smaller scale length to be the bulge. 

Our model consists of two axisymmetric expo- 
nential profiles, so we are unable to fit any non- 
axisymmetric features such as spiral arms, bars, and 
dust lanes. For galaxies in which these features are very 
prominent, the model is inadequate. Additionally, our 
model is inadequate for highly-inclined disks, where the 
semi-minor axis of the galaxy image (i?d x Qd) is compa- 
rable to the disk scale height. At inclinations of ~ 70°, 
the errors in the surface brightness are between 5% and 
10%, for reasonable disk scale heights. If the bulge of 
such a disk galaxy is disk-like (i.e., a pseudo-bulge) but 
is not highly inclined, it is still possible to fit it accu- 
rately. 



4.3 Results 

Panels in Figure [4] shows the observed, model, and 
residual images (from left to right) from B-l-D fits of 
four galaxies (from top to bottom), chosen for illustra- 
tion: SDSSJ013017-hl4 and SDSSJ153433-F41 are disk- 
dominated galaxies, with D/T — 0.99 and 1.00, respec- 
tively; SDSSJ234107— 01 has a more substantial bulge, 
with D/T=0.81; SDSSJ204913-I-00 is an extreme case- 
it has an edge-on disk and a substantial bulge; the disk 
light is obscured by dust, resulting in an extremely low 
reported D/T of 0.54. 

We use both 2.2_Rd and Rso to set the radius at 
which to define the disk rotation velocity amplitude (see 
i6.4|. We note that Rso is less sensitive to any degen- 
eracies in the B-fD fit. The ratio i?8o/(2.2_Rd) clearly 
depends on the relative size of the bulge and disk. Its 
value is 1.4 for a pure exponential disk {D/T = 1), 
and decreases smoothly down to 1.2, 1.1, and 1.0 for 
D/T = 0.8, 0.7, and 0.6, respectively. 

Table|3]lists the best-fit parameters from the i-band 
B-l-D fits of galaxies in the child disk sample, together 



SDSSJ013017 + 14 



with their reduced x"^ values. Figure [s] shows distribu- 
tions of 2.27?d, Rso, qd, and D/T. The 5th, 50th (me- 
dian) and 95th percentile values for 2.27?d are 3.4, 7.4, 
and 12.2 kpc, respectively; for Rso, they are 4.3, 9.8, 
and 16.5 kpc, respectively. 

We find that disk axis ratios from the B-fD fits 
correlate well with the isophotal axis ratios giso from 
SDSS, and are lower by 8 per cent on average. As noted 
earlier, the isophotal axis ratios are affected by seeing 
and bulge light, and both effects tend to increase the in- 
ferred axis ratios. In §5.1[ we describe corrections to qiso 



UN 


% 




SDSSJ153433+41 




SDSSJ234107-01 




SDSSJ204913+00 


UN 







Figure 4. Example two-dimensional bulge-disk decomposi- 
tion fits. Left to right panels; SDSS i-band image, model im- 
age (exponential disk plus exponential bulge, convolved with 
the PSF), and residual image (data— model). The blue ellipse 
on the SDSS image marks the ellipse containing 80 per cent 
of the flux, with semi-major axis Rso and axis ratio equal 
to that of the model disk, q^. Both SDSS and model images 
are logarithmically scaled (with the same stretch), while the 
residual image is linearly scaled. Reading down the figure, 
the maximum (minimum) residual equals 36% (-25%), 18% 
(-11%), 15% (-9%), and 47% (-46%) of the central model 
brightness for each galaxy. 



to estimate disk axis ratios for the parent disk sample 
(to be used for internal extinction corrections). 

Out of the 189 galaxies in the child disk sam- 
ple, 86 (46%), 135 (71%), and 165 (87%) have D/T 
greater than 0.95, 0.9, and 0.8, respectively. We find 
that the majority of galaxies with D/T greater than 
0.95 have ris < 2.3, but there is a large scatter in Us for 
any given D/T. Two galaxies— SDSSJ124545+52 and 
SDSSJ204913-I-00 (shown in bottom panels of Fig.|4|— 
have reported D/T lower than 0.5. Both galaxies have 
edge-on disks {qd < 0.2), which led to an underesti- 
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Figure 5. Distribution of i-band B+D fit parameters- disk 
sizes, 2.2i?(j and ijgo, disls axis ratios, q^, and disk-to-total 
ratios, D/T — for 189 galaxies in the child disk sample. 



mation of the disk light (due to large amounts of dust 
extinction), and a large overestimation of Rso- For both 
galaxies, we do not use the overestimated values for iigo. 
Instead, we assign Rso = 2.3Rd, based on the relation 
between these two disk sizes for B+D model galaxies 
with D/T 0.5. 



5 DERIVATION OF PHOTOMETRIC 
QUANTITIES 

Our aim is to construct photometric estimators of disk 
rotation velocities by calibrating their respective TFRs. 
We begin by describing how we apply internal extinction 
corrections to galaxy absolute magnitudes and colours 
(i5.ll. We do not apply any internal extinction correc- 



tions to the disk scale lengths, but note that their effect 
on the rotation velocity amplitudes used in the TFRs is 
small (see [6.4 for details). 

Next, we define the different photometric quan- 
tities, namely: absolute magnitudes, synthetic magni- 
tudes, stellar masses, and baryonic masses (j ]5.2| -j |575l ). 
We also define and characterize the associated obser- 
vational errors, which are important for estimating the 
intrinsic scatter in the TFRs. In light of future appli- 
cations, we have focused on estimators that are readily 
available for the full parent disk sample (i.e., do not 
require B-fD fits). 



demonstrate that our corrections effectively remove this 
correlation in §10. 1| 

Traditionally, internal extinction corrections are 
applied relative to the face-on orientation, but any in- 
clination (or, equivalently, axis ratio) can be used as a 
reference point. In fact, correcting to a reference incli- 
nation that is typical of the galaxy sample has the ad- 
vantage of minimizing the amount of the correction, as 
well as the corresponding uncertainty. For this reason, 
we use extinction corrections to a reference inclination 
0* (corresponding to an axis ratio q*), which are related 
to face-on corrections by 



(3) 



Here, the superscript ^9 — 9*' denotes corrections to a 
reference inclination 9* , and — 0' denotes corrections 
to face-on orientation. Note that both corrections are 
relative, and do not account for extinction in a face-on 
system. We choose q* = (qd) = 0.40, the mean disk axis 
ratio of the child disk sample]^ 

Following the empirical prescription of |Tully et al.| 



(19981, which is based on a study of 87 spiral galaxies 



in the Ursa Major and Pisces clusters with photometry 
in the Johnson B, Cousins RI, and an infrared K band, 
we calculate the extinction Ax, in the passband A, to be 



-7A(Afr^) X log(gd), 



(4) 



where qd is the disk axis ratio, and M^'~' is the Pet- 
rosian absolute magnitude without internal extinction 
correctionP^ 

The 7a in Eq. [4|are line ar functions of M\ , given in 
Eqs. (3)-(5) of [Tully et all ( |1998| ) for the BRI bands 
(at effective wavelengths 4448, 6581, and 8059A, re- 
spectively). To determine 7a for the SDSS bands ugriz, 
we first use conversion formulae from Table 7 of ISmithl 
|et al.| ( [2002[ ) to estimate BRI absolute magnitudes from 
the SDSS g-band absolute magnitudes, g ~ r and r — i 
colours. Then, we fit a quadratic function to 7b, 7r, and 
7i as a function of A and interpolate (or extrapolate) 7a 
for the SDSS ugriz bands, with effective wavelengths 
3557, 4825, 6261, 7672, and 9097A, respectively (close 
to the BRI passbands). 

For the child disk sample, disk axis ratios qd were 
derived from B-l-D fits (jQ. For the parent disk sample, 
for which we have not performed B-l-D fits, we approx- 
imate the disk axis ratios by applying corrections Cq to 
the SDSS isophotal axis ratios giso 



5.1 Internal extinction corrections 

Disk galaxies are affected by dust obscuration; an in- 
clined disk galaxy appears redder and fainter than if 
it were face-on (e.g., Burstein et al. 19911. Therefore, 



we apply internal extinction corrections to the absolute 
magnitudes and colours of galaxies in both the child and 
parent disk samples. If these corrections are not applied, 
extinction effects can induce a spurious correlation be- 
tween disk inclination and residuals from the TFR. We 



® For comparison, the mean axis ratio of a sample of 
randomly-oriented disks with an intrinsic axis ratio (ratio 
of the vertical and radial scale lengths) of 0.19 is 0.67. Re- 
call that the child disk sample is restricted to low-axis ratio 

galaxies, with qiao < 0.6. 

A more recent work, [Mailer et al. ] | [2009| explored the 
dependence of 7a on both Mn and rig. For rig = 1.0, their 
values of 7a are comparable with those of Tully et al. 1998, 
although the inferred dependence on Mx differs (c.f. their 
Fig. 10). We have not adopted their corrections because they 
require infrared data. 
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Table 3. Best-fit parameters from j-band B+D fits for 189 galaxies in tiie cliild disk sample. 





Galaxy 


-Da(1") 


Id 


D/T 




^80 








(kpc) 






(arcsec) 


(arcsec) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


SDSS 


J001006.61-002609.7 


0.6420 


0.284 


1.00 


3.56 (0.01) 


10.67 


3.91 


SDSS 


J001708. 75-005728.9 


0.3833 


0.575 


1.00 


7.56 (0.04) 


22.62 


2.10 


SDSS 


J002844.82-I- 160058.8 


1.7575 


0.144 


0.83 


5.09 (0.14) 


14.17 


1.23 


SDSS 


J003112.09-002426.4 


0.3943 


0.534 


0.79 


5.38 (0.01) 


14.34 


1.71 


SDSS 


J004916.23+154821.0 


1.5879 


0.568 


0.95 


2.83 (0.12) 


8.28 


1.27 


SDSS 


J004935. 71+010655.2 


0.3570 


0.107 


0.96 


21.93 (0.06) 


64.63 


3.00 


SDSS 


J011750.26+133026.3 


0.6513 


0.430 


0.95 


5.98 (0.08) 


17.53 


1.44 


SDSS 


J012317.00-005421.6 


0.5211 


0.299 


0.81 


5.28 (0.03) 


14.25 


2.75 


SDSS 


J012340.12+004056.4 


0.6654 


0.213 


0.84 


5.55 (0.04) 


15.39 


1.82 


SDSS 


J012438.08-000346.5 


0.5561 


0.221 


0.94 


6.75 (0.03) 


19.76 


1.57 



Notes. — Col. (1): SDSS name of galaxy. Col. (2): Angular diameter distance corresponding to 1" at the redshift of the 
galaxy, in units of kpc. Col. (3): Disk axis ratio, q^. Col. (4): Disk-to-total light ratio, D/T. Col. (5): Disk scale length, R^, 
and its l-cr uncertainty (in parenthesis). Col. (6): Radius containing 80 per cent of the total galaxy light, Rgo- Col. (7): 
Reduced of the B-l-D fit (u is the number of degrees of freedom). 



(jd/giso = Cg(/psF, giso) 



(5) 



where /psp = ^psf/^so, the ratio of the seeing PSF 
FWHM to the galaxy's Petrosian half-light radius. Mo- 
tivated by available data from the child disk sample, we 
adopt the fitting formula 



Cg = 1.21 + mi(/psF - 0.22) + m2{qiso - 0.20), 



(6) 



for galaxies with giso 5^ 0.6. From a least-squares fit to 
the child disk sample dataset, we find mi — 0.22 and 
7712 = —0.32. For galaxies with giso > 0.6, we set Cq — 1. 
The mean value of Cq for the parent disk sample is 1.07. 

By construction, the mean internal extinction cor- 
rections j4®~* are small for the child disk sample: 0.03 
mag for the u, g, and r bands and 0.02 mag for the i 
and z bands. For the i-band, the 5th, 50th, and 95th 
percentile values are —0.22, 0.01, and 0.28 mag, respec- 
tively. The mean internal extinction corrections to face- 
on orientation Al~° are larger: 0.67, 0.60, 0.51, 0.41, 
and 0.31 mag for the ugriz bands, respectively. If we 
were to correct to face-on orientation instead, the main 
effect would be a shift in the derived zero-points of the 
Mx TFRs by an amount given by the relative internal 



extinction corrections (as we show explicitly in [ 9.2 1 



5.2 Absolute magnitudes 

Internal extinction-corrected absolute magnitudes M\, 
as described in 92.11 are calculated as 



:mp,A-51og(DL/10pc)-7^A 



-A^--Ar*.(7) 



Here, mp,A is the apparent SDSS Petrosian magnitude 
in the band A, Dl is the luminosity distance calculated 
from the SDSS redshift (with cosmology fim ~ 0.3, 
f^A = 0.7, and h = 0.7), K\ is the fe-correction to z = 
ca lculated using the kcorre ct product version v4_l_4 



of Blanton & Roweis (|2007l), Af"^ is the correction for 



Galactic extinction based on dust maps oPSchle gel et al.| 
( |l998| ), and ^ is the internal-extinction correction 



to the mean inclination of the child disk sample (Eqs.^ 

&g. - " - 

extinction correction term 



is defined similarly, but without the internal 



We estimate the uncertainty in M\ to be 



5 SVp, 
In 10 cz 



(Mr 



with terms representing the uncertainty in mp^x, -Dl, 
and ^4*"" , respectively. We neglect the smaller un- 
certainties from the redshift, fc-correction, and Galac- 
tic extinction correction measurements. We use un- 
certainties Smp^x derived from the kcorrect rou- 
tines, which have imposed a minimum uncertainty of 
[0.05,0.02,0.02,0.02,0.03] mag in the ugriz bands (to 
account for uncertainties in the absolute calibration for 
each galaxy). For the distance uncertainty, we adopt 
SVpcc = 300kms~^, the typical amplitude of small-scale 
peculiar velocities (Strauss & Willick 19951. The sec- 



ond term has 5th, 50th, and 95th percentile values of 
0.03, 0.07, and 0.12, respectively. For the internal ex- 
tinction correction uncertainty, we follow other authors 



(e.g., P07) and arbitrarily adopt SA^ 



= A". 



/3. 



However, since we correct to the mean inclination (in- 
stead of face-on orientation), our extinction corrections, 
as well as the associated uncertainties, are much smaller. 
The median value of (JA*"" is less than 0.01 mag in all 
bands. For the i band, the 5th and 95th percentile values 
are —0.07 and 0.09 mag, respectively. 

The uncertainty in the absolute magnitude is domi- 
nated by the uncertainty due to peculiar velocities. The 
median uncertainty in Mx is 0.09 mag for the g, r, and i 
bands, and 0.14 mag for the z band. For the i-band, the 
5th and 95th percentile values are 0.04 and 0.15 mag, 
respectively. 
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5.3 Synthetic magnitudes 

We define a synthetic magnitude to be a linear combi- 
nation of the absolute magnitude in some band A and a 
chosen colour, 



A^syn(cf; A, colour) = Mx + a X (colour), 



(8) 



where a is some coefficient that we are free to choose. 
Using this quantity generalizes the method of applying 
colour-based Af,/L corrections to luminosities to esti- 



mate stellar masses (Bell et al. 2003; see S 5.4.1 1 



In this work, we aim to define optimal photometric 
estimators of rotation velocity. Therefore, we choose the 
coefficient a to be the value that minimizes the scatter 
in the inverse TFR, in which Msyn is the independent 
variable, and rotation velocity is the dependent variable 



(c.f. Eq. 26 in f 7.11 



Internal extinction-corrected galaxy colours, as de- 
scribed in 32.11 are calculated as 



^^modcl,r 



model, 3 



A. 



o-e* ^ 



(9) 



where mmodei.A are apparent SDSS model magnitudes, 
-^^modoi,A are fc-corrections to z = for the model magni- 
tudes, A^"^ are corrections for Galactic extinction, and 
are corrections for internal extinction (these last 
three quantities are calculated as in Eq. [?]). {g — r)^'^ 
is defined similarly, but without the internal extinction 
correction terms. 

We estimate observational errors in the g — r colour 
to be 



[5{g - r)]' 



model, r J 



o-e* 



(10) 



Other colours are defined similarly. The median uncer- 
tainty in colours that do not involve the u band range 
from 0.03 to 0.05 mag. Colours involving the u band 
have slightly higher median uncertainties, ranging from 
0.06 to 0.08 mag. 



5.4 Stellar masses 

We consider two kinds of stellar mass estimates: Bell 
et al. (2003) and MPA/JHU stellar masses, defined in 
§5.4. 1| and j ]5.4.2| respectively. We convert both esti- 
mates to correspond to the same IMF normalization. 
We adopt the |Kroupa| ( |2002[ ) IMF, which adequately 
represents direct observational estimates of the IMF for 
late-type galaxies. (The normalization of this IMF is 
0.3 dex lower than that for a Salpeter IMF with a lower 
mass cut-off of 0.1 Mq.) We compare the two different 
stellar mass estimates in |5.4.3[ 



5.4-1 Bell et al. (2003) stellar masses 

Our primary stellar mass estimate, denoted by M« or 
A^,,Bcii, are calculated following the prescription of|Bell| 



et al. (20031. First, we calculate i-band luminosities 

log(L./L0) = -O.4(Mf°-M0,i-hl.lz), (11) 

where M©,; — 4.56 mag, and l.lz is the mean evolution 
correction (in their stellar population models). Then, we 
calculate g — r colour-based stellar mass-to-light ratios 

log(A/./LO = af °" + bf°" x{g- rf" ~ 0.15, (12) 

with af°" = -0.22 and 6f'=" = 0.86. The final 
term, —0.15, converts the normalization from the "diet 
Salpeter" IMF adopted by Bell et al. (2003) to a Kroupa 
IMF Finally, we calculate stellar masses 

log(M./MQ) = log(L,:/L0)+log(M./LO- (13) 

We have used absolute magnitudes and colours 
uncorrected for internal extinction because IBell et aLl 



(20031 did not apply any such corrections. Moreover, as 
noted in Bell & dc Jong (2001), dust affects the stellar 



mass estimates in two different ways that cancel out to 
first order: dimming of the galaxy light leads to a lower 
luminosity and therefore a lower stellar mass estimate, 
but dust reddening leads to a higher stellar mass-to- 
light ratio, and consequently, a higher stellar mass es- 
timate. Because of this partial cancellation, corrections 
for the fact that galaxies in the child disk sample have 
smaller axis ratios compared to the full galaxy popu- 
lation (used by Bell et al. 2003), and therefore larger 
internal extinctions, on average, are tiny, and are not 
applied for simplicity. 

To estimate the statistical uncertainty in A'h, we 
propagate errors from Mi and g — r, as given by Eqs. |8] 
and |10| with the internal extinction correction terms 
excluded. For the child disk sample, the mean statistical 
uncertainty in M«,Beii (at fixed Kroupa IMF) is 0.041 
dex. 



5.4.2 MPA/JHU stellar masses 

The second set of photometric stellar mass estimates, 
denoted by Mt^MPA, are based on fits to the ugriz 
spectral energy distribution (SED) from the MPA/JHU 
groupl^We define M^^mpa to be the reported median of 
the probability distribution function (PDF), with 0.05 
dex subtracted, to convert the normalization from their 
adopted Chabrier IMF to a Kroupa IMF. We define the 
uncertainty in M«,mpa to be half of the difference be- 
tween the 84th and 16th percentiles of the PDF. For 
the child disk sample, the mean statistical uncertainty 
in M*,MPA (at fixed Kroupa IMF) is 0.089 dex. 



The normalization of the "diet Salpeter" IMF is 0.15 dex 
lower than that for a Salpeter IMF with a lower mass cut-ofT 
of 0.1 Mq, and therefore 0.15 dex higher than for a Kroupa 
IMF. 

http:/ /www.mpa-garching.mpg.dc/SDSS/DR7/Data/ 
stellarmass.html. See also Aihara et al.|(i2011|). 
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given by Eq. [13] We adopt u — r colour-based gas-to- 



9.0 9.5 10.0 10.5 11.0 

logM,,Bell(M0) 



Figure 6. Ratio of the two different stellar mass estimates 
A/, Bell S'lid ^t.MPA for the parent and child disk samples 
(yellow dots and red circles, respectively). Error bars show 
l-cr statistical uncertainties in M^, Bell- Dashed and dotted 
curves show l-cr and 2-a bounds for the parent disk sample, 
respectively. 



5.4-3 Comparison of Adt, Ball and Af*,MPA 

Both the Bell et al. (2003) and MPA/JHU steUar masses 
are based on SED fits to SDSS photometry, but they em- 
ploy different methods, including different sets of model 
spectra (the latter allows for a wider possible range in 
parameters) , different stellar population synthesis codes 
( Bruzual fc Charlot|20"03| 



Fioc & Rocca-Volmerange 



19971, and different treatments for extinction (the for- 



mer does not apply internal extinction corrections, while 
the latter does). 

Figure|6]compares the two kinds of stellar mass esti- 
mates for the parent and child disk samples (yellow dots 
and red circles, respectively). There is a clear system- 
atic difference between the two, which increases toward 
lower stellar masses. It is understandable why the low- 
est masses are the most problematic; low-mass galaxies 
have more dust and complicated star formation histo- 
ries. It is somewhat reassuring that the two estimates 
approach agreement at the high mass end, but it is un- 
clear which is more correct for lower masses. 

Given these uncertainties in our current knowledge, 
we will present our main results using both kinds of 
stellar masses. This will give the reader a sense of the 
range of allowed scenarios (at fixed IMF), and explicitly 
show which of our conclusions are robust to systematics 
associated with the stellar mass estimates. 



stellar mass ratio estimates from Kannappan ( 2004 1 
log(Mga./M.,Bdi) = 1.46 - 1.06(u - rf° + 0.15, (14) 



which was derived from 346 galaxies in SDSS DR2 
( Abazajian|2004 1 with infrared data from 2MASS (Jar 



rett et al. 



2000 \ and HI data from HyperLeda ( Paturel 
et al.|2003 1 . We have used colours uncorrected for inter- 
nal extinction because Kannappan ( 2004 1 did not apply 



any such corrections. The addition of 0.15 dex converts 
the normalization of M*,bc11 from their adopted "diet 
Salpeter" IMF to a Kroupa IMF. 

To calculate statistical uncertainties in Mbar, we 
propagate errors from M^. and (u — r)'^*^. For the child 
disk sample, the mean statistical uncertainty in Mbar 
(at fixed Kroupa IMF) is 0.043 dex. 

The above relation, Eq. |14[ has a large scatter, 0.42 
dex, which seems to be mostly physicalFjWe note that 
there are substantial systematic uncertainties in this re- 
lation as well. However, since we are primarily interested 
in using Mbar as a photometric estimator of disk rota- 
tion velocity, its precision as an estimate of the true 
baryonic mass of a galaxy is of secondary importance. 
Regardless, it will be useful to know whether Mbar yields 
a tighter ITER than say, M* or Msyn. 

Figure [7| shows gas-to-stellar mass ratios for the 
parent and child disk samples (yellow dots and red cir- 
cles, respectively). This ratio varies from ~ 30 per cent 
at the highest stellar masses, to a factor of ~ 3 at the 
lowest masses. For the child disk sample, the 5th, me- 
dian, and 95th-percentile values are 0.13, 0.60, and 2.29, 
respectively. 



6 DERIVATION OF KINEMATIC 
QUANTITIES 

We begin this section by describing our long-slit spec- 
troscopy observations and reductions of the raw data 



f 6.1 1. Then, we describe the steps in our rotation curve 



analysis pipeline. First, we extract rotation curves from 



the two-dimensional long-slit spectra (S6.2l. Then, we 



perform fits to these rotation curves (S6.3I. Finally, we 



define the disk rotation velocities that are used in the 



TFRs {^A\. 

In the final subsection, we investigate the effect of 
systematics on the determination of these rotation ve- 
locities. First, we investigate the effect of slit misalign- 
ments (i 6.5.1 1. Next, we compare the rotation curves 



extracted from multiple observations of the same galaxy, 
for the five galaxies observed by both us and P07 



S 6.5.2 1. Finally, we compare the rotation velocities de- 



rived from our analysis pipeline with those reported by 



P07 (i 6.5.3 1 



5.5 Baryonic masses 

We calculate baryonic mass estimates as a sum of stellar 
and gas masses: Mbar = M»,Beii -I- Mgas, with A/»,Beii 



|Kannappan| |2004| obtained a somewhat tighter relation 
using u — K colour, with a scatter of 0.37 dex. Here, we have 
adopted the u — r colour-based relation so as not to rely on 
the availability of K band data. 
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Figure 7. StcUar-to-gas mass ratios derived from m — r 
colours from Kannappan ( 2004 ) (Eq. |14[ | vs. stellar masses 
Bell for the parent and child disk samples (yellow dots and 
red circles, respectively). Error bars show l-cr uncertainties. 
Solid, dashed and dotted curves show the median relation 
and 1- and 2-ct limits for the parent disk sample. 



6.1 Long-slit spectroscopy observations and 
data reduction 

We carried out long-slit spectroscopy observations witli 
tlie 3.5 m telescope at Apaclie Point Observatory (APO) 
using the Dual Imaging Spectrograpli (DIS) witli a 1200 
lines mm~^ grating and 1.5" slit. Tlie spectral range 
covered is 6433-7615 A (centered at 7024 A). Tiie lin- 
ear dispersion is 0.58 A pixel"^, and tiie instrumental 
resolution is approximately 2 pixels or 50kms^^ (at the 
typical observed wavelength of the Ha line). 

We observed a total of 95 galaxies over 25 half- 
nights between March 2009 and June 2010. In most 
cases, we took three 1200 s exposures, except for a few 
bright, nearby galaxies for which we took two 1200 s 
exposures. All targets were observed with the slit po- 
sition angle (PA) set on the major axis of the galaxy, 
determined from SDSS PHOTO exponential fits to the 



r-band images in most cases (see S 6.5.1 for details). Ta- 
ble [4] lists the galaxy names, redshifts, slit PA's, expo- 
sure times, and observation dates. 

To test systematics, we observed two galaxies with 
the slit PA set on the minor axis, and another two galax- 
ies with the slit PA set at ±10 degrees off the major 



axis (see i 6.5.1 1. In addition, five galaxies from the P07 
sample were re-observed with DIS (see [ 6.5.2 l; for these 
galaxies, we use the DIS observations in the final anal- 
ysis. 

Flat-fielding, bias-subtraction, wavelength calibra- 
tion, linearization, and flux calibration were performed 



using standard IRAI^ routines. Wavelength calibra- 
tion and linearization used HeNeAr arc frames. For flux 
cafibration, Feige 34, Feige 56, BD26-I-2606, and HD 
192281 were used as standard stars (selected to match 
the airmass and time of observation). The use of mul- 
tiple exposures for each observation allows cosmic ray 
removal. We combine frames using the IRAF routine 
imcombine with option reject=crreject. For sky sub- 
traction, we took the average sky spectrum over 21 pix- 
els (12 A) on either side of the galaxy's Ha emission. 
Then, from each area of the CCD on either side of the 
galaxy's continuum center, we subtracted the average 
sky spectrum measured from that side. 



6.2 Rotation curve extraction 

We begin by defining the radial bins for the extracted 
rotation curve. First, we choose the central row to be 
at the peak of the galaxy continuum flux (defined to be 
the total flux over rest wavelengths 6590-6630A). Then, 
we bin together rows on either side of the central row 
with increasing bin sizes of 1, 1, 2, 2, 4, 6, 8, 12, 20, 
32, 56, 80, and 100 pixels, for a total of 27 radial bins 
(1 pixel — O.577A). We take the mean over the rows 
in each bin to get binned spectra. With this choice of 
binning, the S/N of the binned spectra remains roughly 
constant with radius (since the bins cover increasingly 
larger areas as the galaxy flux falls off with radius). 

For each of the binned spectra, we fit a set 
of three Gaussians plus a linear continuum to the 
Hq-|-[NII]6548,6583 emission lines, over the wavelength 
range 6520-6610A. We fix the ratio of [NII]6548 to 
[NIIJ6583 flux to the theoretical value 1:3 and require all 
line widths to be equal. We find the best-fit Ha line cen- 
ter Ahq, line width aua, Ha line flux, and [Nil] 6548 line 
flux using a Levenberg-Marquardt least-squares mini- 
mization routine in IDL called MPFIt[^ We perform 
a set of 500 Monte-Carlo realizations (varying the flux 
at each wavelength according to the error in that fiux 
measurement) to estimate errors in the fit parameters. 

For each bin, we also get a flux-weighted average 
spatial position R' , using as weights the total flux over 
the wavelength range Xua ± Saua for each row of that 
bin. Finally, we convert the rest-frame Ha line centers 
and errors into circular velocities to obtain the observed 
rotation curve VobsiRi), with uncertainties av.i- 

As desired, av,i is roughly constant with radius. 
The 5th, 50th, and 95th percentile values for {o'v,i) (av- 
eraged over all radii) are 0.86, 2.43, and 9.60 kms~^, 
respectively. Splitting each rotation curve at its turn- 



over radius -Rto (defined in Eq. 15 below), the {(Tv,i), 
averaged over the inner and outer regions, have median 
values of 2.96 and 2.10 kms~^, respectively (the 5th 
and 95th percentile values are similar to the full radial 
average as well). 



http://iraf.noao.edu 
http://purl.com/net/mpfit 
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Table 4. Observations with DIS at APO 3.5m (95 galaxies). 



Galaxy name 


z 


PA 


Exp. time 


Obs. Date 






(deg) 


(sec) 




(1) 


(2) 


(3) 


(4) 


(5) 


SDSS J002844.82+160058.8 


0.0947 


35.79 


1200 X 3 


2009-09-14 


SDSS J004916. 23+154821.0 


0.0846 


-73.04 


1200 X 3 


2009-09-14 


SDSS J011750. 26+133026.3 


0.0326 


-66.45 


1200 X 2 


2009-09-15 


SDSS J013017.16+143918.5 


0.0239 


-23.95 


1200 X 3 


2009-09-23 


SDSS J020056. 00+133116.6 


0.0312 


41.92 


1200 X 3 


2009-09-23 


SDSS J020133. 02+133126.4 


0.0978 


48.43 


1200 X 2 


2009-09-14 


SDSS J020526. 65+131938.4 


0.0251 


78.57 


1200 X 2 


2009-09-20 


SDSS J020540. 31-004141.4 


0.0424 


-30.32 


1200 X 2 


2009-09-15 


SDSS J020819. 43+134944.6 


0.0798 


-22.48 


1200 X 3 


2009-09-23 


SDSS J020923. 14+125029.8 


0.0604 


64.47 


1200 X 2 


2009-09-15 



Notes. — Col. (1): SDSS name of galaxy. Col. (2): Redsliift of galaxy. Col. (3); Major axis position angle in degrees East of 
North. Col. (4): Duration of spectroscopic exposure (sec). Col. (5): UT date of observation. 



6.3 Rotation curve fitting 

Studies have shown that most observed disk galaxy ro- 
tation curves can be modeled by an arctangent model 
( Courteau|1997| |: 

d(i?') = Vo + -K.obsarctan ( ^ . (15) 



Rto 



This model has four free parameters: the systemic veloc- 
ity Vo, the asymptotic circular velocity K,obs, the spa- 
tial center Rq, and the turn-over radius Rto, at which 
the rotation curve starts to flatten out. We denote the 
radius defined from the centre of the rotation curve as 
R = R' — Ro- Note that all the velocities above are 
not corrected for inclination. 21 Using the Levenberg- 
Marquardt routine MPFIT in IDL, we find the best-fit 
parameters by minimizing x^, defined as 



2 _ [VmodiR'i) — VohsiR'i)]^ 



(16) 



where the sum is over radial bins with sufficiently 
good S/N spectra {S/N > 5 in the Ha line), and 
the effective uncertainty in Vohs{R'i) is defined to be 
foff,! = {cTyi + (Tv_add)^^^- Here, av,i is the formal un- 
certainty from the Gaussian line fits (described in [ 6.2 1, 



and (Tv,add is an additional model uncertainty that ac- 
counts for non-circular motions in the galaxy. We set 
o"v,add = lOkms^^, following P07, but note that the 
fits are not very sensitive to this choice. Averaging over 
all radii, {(Jeff ,i) has 5th, 50th, and 95th percentile values 
of 10.05, 10.37, and 14.89 kms"\ respectively. 

We define the total uncertainty in a fit parameter 
P (e.g., Vc.oha, Rto, etc.) as 



1/2 



(17) 



where the fit uncertainty am is the formal error from 
the minimization fit and the systematic uncertainty 
(Jsya is given by 



where P'"^' is the result of repeating the fit with a radius- 
weighted scheme that gives greater weight to points at 
larger radii, i.e., replacing the denominator in Eq. [16] 



with = {l/\R'\)a! 



By adopting this definition of total uncertainties, to 
be used in our fits to the TFRs, we automatically down- 
weight the contribution of those galaxies for which the 
flat part of the rotation curve is not well-constrained 
(e.g., irregulars, galaxies in mergers, and other pecu- 
liar cases). A large fit uncertainty indicates that Vc,obs 
is not well-constrained because the S/N is low, or be- 
cause the observed rotation curve does not show a turn- 
over (i.e., it is still rising at the outermost point). On 
the other hand, a large systematic uncertainty indicates 
that Vc,obs is not well-defined because the arctangent 
model is not a good model for the observed rotation 
curve. 

For illustration, Figure |8] shows the unweighted 
and radius-weighted fits (red solid and blue dashed 
curves, respectively) for two galaxies. The top panel 
shows the observed rotation curve of SDSSJ124259+42, 
a case wherein the two fits coincide, resulting in a small 
systematic uncertainty in 14, obs- The bottom panel 
shows the peculiar rotation curve of SDSSJ170952+35, 
one of the few cases wherein the two fits disagree 
considerably, resulting in a large systematic uncer- 



sys 



(K:,obs)/K:,o' 



bs 



(18) 



tainty. In this extreme case, 

|(-157 kms-^)/(135 kms'^)] = 1.16. By compar- 
ison, the 5th, 50th, and 95th percentile values of 
|crsys(K,obs)/K,obs| for the cliild disk sample are 0.0035, 
0.033, and 0.14, respectively. 

Table [5] lists parameters from the arctangent model 
fits for the 189 galaxies in the child disk sample. We 
list the turn-over radii -Rto and inclination-corrected 
asymptotic circular velocities K and rotation velocity 
amplitudes V2.2 and Vso (defined in S6.4I. Rto is listed 
with total 1-a uncertainties, while all rotation velocities 
are listed with both the total and formal fit 1-a uncer- 
tainties (the latter are inside parentheses). The table 
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Galaxy name 


sin 8 


Rto 


Vc 


V2.2 


Vao 


y? 1^ 






(arcsec) 


(kms^"*") 


(kms-i) 


(kms^"*") 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


SDSS J001006.61-002609.7 


0.9765 


2.66(0.47) 


158.2 ± 10.8(5.7) 


125.2 ± 5.7(1.0) 


133.6 ±6.9(2.1) 


2.7 


SDSS J001708. 75-005728.9 


0.8336 


2.81 (0.98) 


137.8 ± 16.9(12.8) 


123.1 ± 13.5(6.9) 


126.9 ± 14.1(8.4) 


2.8 


SDSS J002844.82+160058.8 


1.0000 


1.45(0.68) 


144.1 ± 16.7(15.2) 


132.3 ±9.9(8.2) 


134.7 ± 11.3(9.6) 


1.5 


SDSS J003112.09-002426.4 


0.8611 


0.76(0.22) 


153.9 ±9.3(7.3) 


147.7 ±8.6(5.5) 


148.8 ±8.7(5.8) 


0.8 


SDSS J004916.23+154821.0 


0.8386 


1.67(0.34) 


150.9 ± 12.2(1.9) 


125.7 ± 13.4(0.7) 


131.8 ± 12.9(1.0) 


1.7 


SDSS J004935. 71+010655.2 


1.0000 


17.70(1.52) 


267.2 ± 11.0(1.3) 


207.4 ± 14.4(0.4) 


221.8 ± 13.9(0.6) 


17.7 


SDSS J011750.26+133026.3 


0.9194 


4.06(0.92) 


108.8 ± 11.5(2.4) 


88.1 ± 7.9(1.6) 


93.1 ± 9.0(1.8) 


4.1 


SDSS J012317.00-005421.6 


0.9718 


2.23(0.57) 


197.7 ± 13.1(11.4) 


173.8 ± 7.2(4.8) 


178.2 ±8.3(6.0) 


2.2 


SDSS J012340.12+004056.4 


0.9952 


1.76(0.24) 


206.6 ±8.6(6.1) 


187.8 ±6.1(3.6) 


191.6 ±6.6(4.1) 


1.8 


SDSS J012438.08-000346.5 


0.9933 


4.18(0.42) 


196.7 ± 7.3(2.3) 


162.3 ±9.4(0.6) 


170.6 ±9.1(1.0) 


4.2 



Notes. — Rotation velocities (cols. 4-6) are listed with their l-cr total (after the ± symbol) and formal fit uncertainties (in 
parentheses). Col. (1): SDSS name of galaxy. Col. (2): Inclination correction. Col. (3): Turn-over radius, with l-cr total 
uncertainty. Col. (4): Inclination-corrected asymptotic circular velocity. Col. (5): Inclination-corrected velocity amplitude at 
2.2/?d. Col. (6): Inclination-corrected velocity amplitude at -RgQ. Col. (7): Reduced of the fit. 
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Figure 8. Observed rotation curves for SDSSJ124259+42 
and SDSSJ170952+35 (top and bottom panels, respectively), 
shown with the best-fit arctangent models from both un- 
weighted and radius-weighted schemes (solid red and blue 
dashed curves, respectively). Vertical lines mark character- 
istic radii, i? = 0, 2.2i?(j, and ijgo (dotted, dot-dashed, and 
dashed lines, respectively). 



also lists inclination corrections sinS and reduced 
values of the fits. 

The differences between the fit x^ values are not 
being driven by statistical (photon) noise, so their actual 
values do not carry the expected meaning. Instead, they 
are driven by systematic deviations from the arctangent 
model, due to asymmetry or "hooks" in the observed 
rotation curve. So although the typical value of 'x^ jv is 
0.25, a value larger than unity does not imply that the 
observed rotation curve is not globally well-described by 
the arctangent model. 

Figures [9| jll| show folded, normalized, and scaled 
rotation curves and best-fit arctangent models for the 



189 galaxies in the child disk sample. The rotation 
curves have been folded at the best-fit spatial centre 
i?o, normalized by the asymptotic circular velocity 14, 
and scaled by the disk scale length R^. Panels are la- 
belled with a shortened galaxy name and -Rto IRa - In all 
but two cases (1 per cent of the sample), the observed 
rotation curve can be adequately modeled on average 
by the arctangent model, albeit many galaxies show 
hooks. The only two exceptions are SDSSJ170952-f 35 
(also shown in Fig. [8| and SDSSJ141026-00, for which 
the observed rotation curve turns over toward lower ve- 



locities instead of flattening out (see Fig. 10 1. The flts 
to these two galaxies have extremely poor x values 
of 2.5 and 21.6, respectively. 

For the child disk sample, the median value of the 
total \-a uncertainty in log 14 is 0.030 dex. As expected, 
this uncertainty decreases toward higher stellar masses 
(i.e., higher SjN data). The median values are 0.045, 
0.029, and 0.020 dex for the bottom, intermediate, and 
top 1/3 bins in M*, respectively. If the systematic un- 
certainty is not included, the median value decreases to 
0.024 dex for the fuU sample, and to 0.037, 0.024, and 
0.015 dex, for the three M» bins, respectively. 

Figure [12] shows the relation between stellar mass 
M« and the ratio of turn-over radius to disk scale 
length, RTo/Rd- We confirm that the most luminous 
galaxies tend to have very steeply-rising rotation curves 
(e.g., [Persic et aL][T996} [Catinella erar]|2006| . At the 
other extreme, we find two outliers with RTo/Rd > 5 
(SDSSJ080046+35 and SDSSJ120626+24). The rota- 
tion curves of both galaxies are observed out to their 
i?80, but are still rising at the outermost point[^ Con- 
sequently, the asymptotic circular velocity is not well- 
constrained, resulting in extremely large values of a{Vc) 

SDSS J 120626+24 appears to be a central galaxy of a 
group, which explains its rising rotation curve beyond its 
optical radius. 
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Figure 9. Normalized, folded, and scaled rotation curves for 63 (out of the 189) galaxies in the child disk sample, sorted by 
galaxy name. Velocities are scaled by the best-fit asymptotic circular velocity Vc; radii are scaled by the disk scale length ijj 
and shown in a logarithmic scale from 0.05— 5i?d- Filled and open circles show opposite arms of the rotation curve, folded at 
the best-fit spatial centre Rq. The solid curve shows the best-fit arctangent model. Dashed and dotted lines mark Rto E^nd 
i?di respectively, and the vertical bar marks Rgo, the radius containing 80 per cent of the i-band galaxy light. Note that we 
define rotation velocity amplitudes at optical radii 2.2Hd and Rgo- Panels are labelled with the galaxy name and the ratio of 
the turn-over radius to the disk scale length, R-po/^d- The rest of the galaxies are shown in Figs. [l0||ll| 
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Figure 10. Same as Fig.^ Normalized, folded, and scaled rotation curves for 63 galaxies in the child disk sample (2 of 3). 
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Figure 11. Same as Fig.^ Normalized, folded, and scaled rotation curves for 63 galaxies in the child disk sample (3 of 3). 
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Figure 12. Relation between the ratio of turn-over radius 
to disk scale length RTo/^d ^.nd stellar mass Af, bcII ior 
189 galaxies in the child disk sample. Two outlier galaxies — 
SDSSJ080046+35 and SDSSJ120626+24— have RTo/Rd > 
5 and lie beyond this plot. 



of 960 kms"""^ and 605 kms"'^, respectively, as well as 
largely overestimated values of Vc, equal to 398 kms^^ 
and 457 kms~^, respectively. 



symmetry 1^ 

Vc = K.obs/sinS, 
V2.2 = V2.2,obs/sin6l 
V80 = 14o,obs/sin6l = 



-- V{R = 2.2Rd)/sme, 
V{R^ R8o)/sme. 



Here, the inclination corrections are given by 



1 



1 



1/2 



(19) 
(20) 
(21) 



(22) 



where is the intrinsic axis ratio, or the ratio of the 
vertical and radial scale lengths of the disk. Following 
P07, we adopt a single value for = 0.19 ( Haynes fc 
GiovaneTIIl|1984 1, although it is known to vary slightly 
from disk to diskF^ 

Figure [13] shows cumulative histograms in the ra- 
tios -RoptZ-Riast and Kot/Viast, where _Riast is the radius 
of the outermost point in the observed rotation curve, 
and Viast is the circular velocity at that radius based on 
the best-fit arctangent model. We find that most of the 
observed rotation curves extend close to or beyond Rso, 
so a substantial extrapolation of the rotation curve is 
rarely required. Moreover, Rso is far enough out to be 
close to Vc in most cases. This is also evident in the pan- 
els of Figs. [9[ [11| where Rso is marked by solid vertical 
bars. These results support the empirical logic of adopt- 
ing the rotation velocity amplitude Vso (as in P07). 

Moreover, Vgo is relatively insensitive to changes 
in the value of Rso because it samples the flat portion 
of the rotation curve in most cases. For example, even 
though Rso varies from 1.2-1.4 x(2.2i?d), Vso is within 
10% of V2.2 for most galaxies, and is only 5% larger on 
average. Therefore, the uncertainty in Rso or Rd, in- 
cluding, e.g., the effect of internal extinction]^ has a 



6.4 Rotation velocity amplitudes 

We now define inclination-corrected rotation velocity 
amplitudes that we will use in the TFRs. In addition 
to the asymptotic circular velocity 14, we define rota- 
tion velocities evaluated at some suitably chosen optical 
radius -Ropt. A common choice for -Ropt is 2.2 times the 
disk scale length 2.27?d ( Courteau fc Rix|1999 l, chosen 
because this is the radius at which the rotation curve 
of a pure self-gravitating exponential disk would peak 
( Freeman|1970 1. We also adopt Rso, following P07, the 
radius containing 80% of the i-band flux. As noted in 
S4.3 Rso/Rd = 3.03 for a pure exponential disk, but 
is smaller for galaxies with significant bulges. Both Rd 
and -Rso have been corrected for the effect of seeing, but 
compared to Rd, Rso has the advantage of being less 
sensitive to the degeneracies of the B-fD fits. Moreover, 
Rso is closer to the peak of the total (disk-|-DM halo) ro- 
tation curve, which is located at ~3-Rd for typical haloes 
of disk galaxies ( |Mo et al.|1998[ ). 

We follow convention and correct the disk rotation 
velocities to the edge-on orientation (assuming circular 



Without assuming circular symmetry, the rotation veloc- 
ity amplitudes would depend not only on inclination, but 
also on the orientation of the galaxy with respect to the line 
of sight. When we look down the long axis of gas orbits, the 
observed velocities are larger than when we look down the 
short axis. Deviations from circular symmetry and the vari- 
ation in viewing angles contribute to the observed scatter in 



the TFRs (see {7.3 for further discussion). 

From analyses of the distribution of axis ratios of SDSS 
galaxies , |Ryden| l |2004| found = 0.22 ±0.06 and |Padilla"fcl 



Strauss 



2008 1 found g, 
authors 



0.21 ± 0.02. 

have corrected the observed disk scale 



lengths for internal extinction. We do not apply such cor- 
rections because the magnitude and even the sign of the cor- 
rection, which depends on the spatial distribution of dust in 
the disk, is a priori unclear. [Mailer et al.| | |2009[ l found that 
more highly-inclined late-type galaxies have larger half-light 
radii than face-on ones, indicating that attenuation of dust 
is stronger in the center of galaxies. [Courteau et al.] l |2007^ 
applied a correction factor of (1 + 0.4 log(a/f)))~^ , following 
'Giovanelli et akj l [1997[ |. This correction is at most 29 per 
cent, comparable to the difference between Rso and 2.2i?,d. 
The associated systematic bias in Vso is thus comparable to 
the typical measurement uncertainty in Vso, of around 0.03 
dex. 
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Figure 13. Cumulative distributions of i?8o/^?iasti 
2.2i?(j//?iast (thick and thin curves in the upper panel, 
respectively) and Vgo/Mast ^.nd V2.2/Viast (thick and thin 
curves in the lower panel, respectively) for the child disk 
sample. 



relatively small effect on the rotation velocity ampli- 
tudes. 

To calculate the total error in Vso, we propagate 
the formal errors in the best-fit model parameters using 
the full covariance matrix, then include a systematic 
uncertainty similar to that defined for the fit parameters 



(Eq. 18 1, as well as a contribution from the uncertainty 



in the inclination correction 



(23) 



R=R8a 



(24) 



Here, V{R) is the arctangent model function defined in 
Eq. |15| and ai = {Vq, Vc, iio, ^to} are the fit param- 

(t) 

eters with formal fit covariance matrix dj, Vgg is the 
velocity amplitude at Rso for the best-fit model derived 



from a fit using the radius- weighted scheme (c.f. Eq. 18 1, 
and 5{sm6) is based on the formal fit error in qd- &V2.1 
is defined similarly. 

The error in Vgo is largely dominated by the for- 
mal fit error; the contribution from the formal inclina- 
tion correction uncertainty is negligible. The systematic 
uncertainty in Vso is smaller than that in 14, by ap- 
proximately half. For the child disk sample, the 5th, 
50th, and 95th percentile values of log Vgo are 0.012, 
0.027, and 0.10 dex, respectively. The errors in logV2.2 
are similar. 



6.5 Tests of systematics 

We now investigate potential systematic effects on the 
measured rotation curves. In §6.5. 1[ we look at the ef- 
fect of slit misalignments by comparing observations at 
different slit PAs (along the minor axis and ±10° off the 
major axis). In i 6.5.2 we look for any systematic dif- 



ferences between observations performed by P07 with 



CCDS/TWIN and those performed with DIS. In j |6X3 
we compare the best-fit parameters from our rotation 
curve fits with those from P07. 



6.5.1 Off-axis PA observations 

We do not have a direct measurement of the kinematic 
major axis PA, and rely on photometry to determine 
the major axis PA. The SDSS photometric pipeline pro- 
vides three different measurements of the photometric 
major axis PA: (i) from the exponential and deVau- 
coleurs model fits, (ii) from ellipticity measurements 
at the 25 mag arcsec"^ isophote, and (iii) from fiux- 
weighted adaptive second moments. The B-l-D fits per- 
formed in ij4] provide another measurement of the disk 
major axis PA. The model and B-l-D fits account for the 
effect of the seeing, but the other two measurements do 
not. 

If the velocity field is noncircular, the photomet- 
ric major axis may be offset from the kinematic major 
axis. For a potential with an ellipticity of 0.1, the mean 
PA difference is 5° (averaged over all viewing angles), 
roughly in agreement with observed values ( Franx & de 
Zeeuw|[l992 |. In addition, there is some uncertainty in 
the determination of the photometric major axis PA it- 
self. Comparing the different photometric major axis 
PA estimates, we find that except for one outlier — 
SDSSJ093724-I-48— all galaxies in the child disk sam- 
ple have an isophotal PA within ±11°, and a B+D fit 
PA within ±6°, of its exponential fit PA. Taking into 
account both effects, the typical level of PA offsets is 
expected to be < 10°. 

As the default configuration, we observed galaxies 
at the photometric major axis PA determined from the 
r-band exponential model fit performed by the SDSS 
photometric pipeline. This was done for 86 out of the 
95 targets, including SDSSJ093724+48. Seven galax- 
iej^ were observed at the r-band isophotal PA, and 
two galaxiej^ were observed at a PA determined by 
eye, which turn out to be within 3° of the exponential 
fit PA (and within 1° of the isophotal PA). All observed 
slit PAs are listed in Table HI 

For the outUer SDSSJ093724±48, a two-armed spi- 
ral, the ideal slit orientation is unclear (see left panel 
of Fig. 141. We first observed at the exponential fit PA 



20 These galaxies are: SDSSJ110148±44, SDSSJ124021-I-48, 
SDSSJ131159±43, SDSSJ132907±41, SDSSJ133151±55, 
SDSSJ143117±47, and SDSSJ171157+58. 

21 These galaxies are SDSSJ142201+44 
SDSSJ151434-I-39. 



and 
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SDSSJ093724+48 




R (arcsec) 

Figure 14. Left panel: SDSS image of the galaxy SDSS J093724+48, the only case where the exponential fit and isophotal PAs 
are significantly different. We observed this galaxy at two slit positions: the r-band exponential fit PA (E of N) of 15.88° (blue 
rectangular box) and at an eye estimate of the galaxy major axis PA of 60.9° (green rectangular box). The isophotal PA of 
70.9° is also shown for reference (thick red line). In this image, North is up and East is toward the left; the horizontal line on 
the lower-right corner is 5" long; the green box is 1.5" by 40" (the actual slit covers 1.5" by 2', so its length spans the whole 
extent of the galaxy). Right panel: Rotation curves observed at the two slit PAs shown in the left panel. The observation at PA 
= 60.9° yields a best-fit asymptotic circular velocity 3.6 times larger than that for the observation at PA = 15.88° ((red circles 
and solid curve vs. blue stars and dashed curve), a strong indication that this PA is closer to the true major axis of the galaxy; 
therefore, we adopt it in our final analysis. 







r _■ 







Figure 15. SDSS images of four galaxies with off-axis slit PA observations. Here, green rectangular boxes indicate the major 
axis PAs, and blue boxes indicate the off-axis slit PAs. As in Fig. |14| North is up and East is toward the left; the horizontal line 
on the lower-right corner is 5" long; green rectangles are 1.5" by 40". Upper panels: SDSSJ144658-I-12 and SDSSJ154014-I-50 
were observed along the minor axis PA (blue rectangular boxes). Lower panels: SDSSJ020526-I-13 and SDSSJ130140-I-08 were 
observed at ±10° off the major axis PA (blue rectangular boxes). 
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of 15.88° (blue rectangular box in Fig. 14 1, and then at 
an eye estimate for the galaxy major axis PA of 60.9° 
(green rectangular box). The SDSS isophotal PA is 70.9° 
(thick red line), 10 deg off the observed PA. The right 
panel of Fig. |14| shows the rotation curves and best-fit 
arctangent models for the two observations, obtained 
using procedures described in §6.2| and §6.3[ The slit 
orientation with PA= 60.9° clearly yielded a larger ro- 
tation velocity amplitude, with a best-fit asymptotic cir- 
cular velocity (before inclination correction) of K,obs = 
148±16kms~^, compared to 41±7kms~^ for the other 
slit orientation. This orientation is therefore closer to 
the true major axis of the galaxy, and we adopt this 
observation in our final analysis. 

To check for a velocity gradient along the mi- 
nor axis, a clear signature of noncircular motions, we 
observed two of our targets — SDSSJ144658-I-12 and 
SDSSJ154014-I-50 — along the minor axis, i.e., perpen- 
dicular to the observed major axis PA (blue and green 
rectangular boxes in the upper panels of Fig. [15] show 



the minor and major axis PAs, respectively). Figure 16 
compares the minor and major axis rotation curves and 
best-fit arctangent models (red circles/solid curves and 
blue stars/dashed curves, respectively). 

For SDSSJ144658-I-12, the best-fit parameters are 
Vc.obs = 181 ± 17kms-i and Rto = 3.8 ± 0.6", 
with a reduced of l-l- Fitting the minor axis ro- 
tation curve with -Rto fixed at this value, we find 
Vc.obs = 18 ± 29kms~^, consistent with zero (with a 
reduced of 0.25). Performing a similar procedure for 
SDSSJ154014-h50, we find K.obs = 151 ± 6 and 19 ± 
6kms~^, with reduced values of 0.57 and 0.32, for 
the major and minor axis observations, respectively. In 
this case, there is some evidence for rotation in the mi- 
nor axis rotation curve, but only at 1 /8 of the value of 
that for the major axis rotation curve. 

To test the effect of slit misalignments on the 
inferred velocity amplitudes, we observed two of our 
targets— SDSSJ020526+13 and SDSSJ130140-h08— at 
±10° offsets from the observed major axis PAs (blue 
rectangular boxes in the lower panels of Fig. [15] show 
the offset slit PAs and green boxes show the major axis 
PA) . Figure [Tt] compares the observed rotation curves 
for the different slit orientations and the best-fit arc- 
tangent models. For both galaxies, one of the off-axis 
rotation curves has a slightly lower velocity amplitude 
compared to the other two (as expected if the major axis 
PA is well determined). Table [6] lists the best-fit arct- 
angent parameters for these observations, together with 
their formal 1-a fit uncertainties. The rotation velocity 
amplitudes 14o,obs (col. 6) for the different observations 
are consistent to within la, indicating that this quan- 
tity is robust to slit PA offsets of order 10° (the expected 
level of slit misalignments, as noted earlier). 

We conclude that slit misalignments are not ex- 
pected to introduce a significant systematic bias in the 
TFR toward lower rotation velocity amplitudes. We 
note, though, that they are expected to contribute to the 
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Figure 16. Major and minor axis rotation curves for 
SDSSJ144658+12 and SDSSJ154014+50 observed along the 
minor axis, i.e. at 90° off the major axis (blue stars and 
dashed curves), compared to those observed along the major 
axis (red circles and solid curves). 



6.5.2 Comparison with P07 observations 

Since the child disk sample combines long-slit spec- 
troscopy observations taken using different instruments, 
it is important to check that this difference does not in- 
troduce any systematic bias in the derived kinematic 
quantities. For this purpose, we have randomly selected 
five out of the 99 galaxies from the P07 sample that 
satisfy our selection criteria, and re-observed them with 
DIS. 

Figure |18] compares the rotation curves and best- 
fit arctangent models for these galaxies. There is a good 
match between the different observations (blue stars and 
dashed curve for TWIN or CCDS, and red circles and 
solid curve for DIS). For each galaxy, the best-fit ro- 
tation curve from the DIS observation is consistent, to 
within a few a, with the observed rotation velocities 
from the TWIN or CCDS observation, and vice versa. 

Table [7] lists the best-fit arctangent model parame- 
ters, together with their 1-a formal fit uncertainties, for 
the different observations. The rotation velocity ampli- 
tudes Vgo.obs (col. 6) for the different observations are 
consistent to within la for four out of the five galaxies. 
The exception is SDSSJ020540-00, for which the arct- 
angent fits are affected by hooks in the rotation curve 
(evident around 5"). Focusing on the outer arms, it is 
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Table 6. Rotation curve fits for targets with both major axis and off-axis observations (2 galaxies). 



Galaxy name 


Slit PA 






Rto 




VgO.obs 






(deg) 


(km s~ 


h 


(arcsec) 


(arcsec) 


(kms"-*-) 




(1) 


(2) 


(3) 




(4) 


(5) 


(6) 


(7) 


SDSSJ020526-I-13 


78.57 


144.9 ± 


5.0 


2.43 ±0.26 


25.27 


136.0 ±4.8 


0.27 




68.57 


149.6 ± 


5.3 


2.69 ±0.29 




139.5 ± 5.1 


0.36 




88.57 


137.1 ± 


5.4 


2.58 ±0.31 




128.2 ± 5.1 


0.22 


SDSSJ130140+08 


-45.43 


149.8 ± 


6.4 


1.78 ±0.21 


10.49 


133.8 ±6.0 


0.41 




-55.43 


136.6 ± 


6.7 


1.72 ±0.23 




122.4 ±6.3 


0.71 




-35.43 


157.1 ± 


8.1 


2.50 ±0.32 




133.7 ±7.4 


0.26 



Notes. — Col. (1): Galaxy name. Col. (2): Slit position angle in degrees East of North. The first line is the major axis PA; the 
next two lines are offsets by ±10°. Col. (3): Best-fit observed asymptotic circular velocity \4,obs- Col. (4): Best-fit turn-over 
radius -Rto- Col. (5): Radius containing 80 per cent of the i-band light. Col. (6): Velocity amplitude evaluated at -Rgo, before 
inclination correction, Vgo.obs- Col. (7): Reduced of the fit. Cols. (3), (4), & (6) are listed with their formal 1-cr fit errors 
(not including systematic uncertainties). 
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Figure 17. Rotation curves for SDSSJ020526+13 and 
SDSSJ130140+08 observed at ±10° off the major axis (green 
triangles and dot-dashed curves, and blue stars and dashed 
curves, respectively), compared to those observed along the 
major axis (red circles and solid curves). 



clear that the rotation curve data themselves are con- 
sistent with one other. 

In detail, the shapes of the rotation curves dif- 
fer from one another, as reflected in the differences 
in the best-fit Rto- This is especially apparent for 
SDSSJ024459+01 and SDSSJ152137+00, in which the 
differences in Rto are 50 and 40 per cent, respectively. 
This difference can be attributed to differences in see- 



ing conditions. We conclude that there is no evidence for 
any systematic bias between the long-slit spectroscopy 
observations performed by us and P07. 



6.5.3 Comparison with P07 rotation curve fits 

In this section, we compare our best-fit arctangent 
model parameters — Vc,obs and Rto — with those pub- 
lished in P07, for 157 galaxies in the P07 sample (99 of 
which are in our child disk sample) . We use the long-slit 
spectroscopy observations of P07, but perform our own 
analysis of the data to extract and fit rotation curves. 
Comparison of our results with theirs provide an esti- 
mate of the level of systematic differences that can arise 
from different analysis methods. 

Figure [19] shows the differences between P07 mea- 
surements and ours, A(logX) — (logX)^^*"^' — logX, 
where X — {Rto, Vc.obs, Vso, sin^} (top to bottom pan- 
els). Nineteen out of the 157 galaxies have |AVc,obs| > 
2o-(\4,obs) (marked by red circles, in all panels), and 
for almost all of them, AVc.obs < 0. We surmise that 
the reason for this is that our extracted rotation curves 
tend to extend out farther than those of P07, which led 
to lower (i.e., less overestimated) K.obs- As expected, 
we find that the differences in Vgo are smaller than 
those in Vcobs, explicitly illustrating the relative robust- 
ness (and therefore, advantage) of this choice of veloc- 
ity amplitude. Finally, we note that differences in the 
estimated disk axis ratios led to differences in Vso as 
large as ~0.05 dex (see bottom panel). Both we and 
P07 used two-dimensional bulge-disk decompositions to 
determine disk axis ratios, but our fits were performed 
on images from a more recent SDSS reduction. We use 
our own fitting routine (described in ^ , while P07 u sed 
the publicly-available code GALFIT |Peng et al.|2002j ). 

Later in §9.5| we compare the TFRs derived from 
our and P07's measurements. We find that the M—Vgo 
relations are consistent within the reported 1-cr uncer- 
tainties, indicating that our results are robust to sys- 
tematic differences between the analysis pipelines. 
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Figure 18. Comparison of rotation curves obtained from DIS (red circles) and P07 TWIN and CCDS observations (blue stars). 
Best-fit arctangent models are shown in red solid and blue dashed curves, respectively. Dot-dashed vertical lines mark Rgo, the 
radius enclosing 80 per cent of the i-band light, at which the rotation velocity amplitude Vgo is defined. 



7 DERIVATION OF TFRS 

In this section, we describe our derivation of the scal- 
ing relation between the disk rotation velocity Kot and 
some photometric quantity O. The difTerent photomet- 
ric quantities we consider — absolute magnitudes M\, 
synthetic magnitudes Msyn, stellar masses M*^Beii and 
M, MP A, and baryonic masses Mbar — were defined in 

We aim to identify the "optimal" photometric esti- 
mator of disk rotation velocity, calibrate it, and charac- 
terize its scatter, using the child disk sample (189 galax- 
ies). The photometric estimates can then be applied to 



the full (spectroscopic) parent disk sample (~ 170 000 
galaxies), and potentially, to an even larger photometric 
galaxy sample. 

We describe the modeling and fitting of the TFRs 
in §7.1| and |7.2[ respectively. In j ]7.3| we describe how 
the fits should be interpreted. Results of the fits will be 
presented in the next two sections, i|8]and ij9] 

7.1 Modelling the TFRs 

Since we are mainly interested in calibrating pho- 
tometric estimators of disk rotation velocities, we 
will focus on fits to the so-called inverse TFRs (IT- 
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Table 7. Rotation curve fits for targets with observations from both P07 (CCDS/TWIN) and DIS (5 galaxies). 



Galaxy name 


Instru. 


^^,obs 


Rto 










(deg) 


(km s~-^) 


(arcsec) 


(arcsec) 


(km s"'^) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


SDSSJ020540-00 


DIS 


216.2 ± 3.6 


0.80 ± 0.06 


29.13 


212.4 ± 3.5 


8.33 




TWIN 


201.7 ±3.5 


0.85 ± 0.07 




198.0 ±3.5 


10.11 


SDSSJ024459+01 


DIS 


221.9 ± 10.2 


4.93 ± 0.49 


19.35 


186.6 ±9.2 


0.36 




TWIN 


206.3 ±8.8 


3.20 ± 0.32 




184.8 ±8.1 


0.53 


SDSSJ080927+43 


DIS 


170.1 ± 5.6 


1.77 ± 0.21 


15.66 


157.9 ±5.2 


1.31 




CCDS 


184.4 ±7.2 


2.42 ± 0.36 




166.4 ±6.0 


0.38 


SDSSJ084408+50 


DIS 


125.3 ±9.5 


4.83 ± 0.83 


21.12 


107.4 ± 11.1 


0.36 




CCDS 


121.9 ± 10.4 


5.72 ± 1.13 




101.4 ± 13.3 


0.26 


SDSSJ152137+00 


DIS 


155.6 ±8.1 


3.29 ± 0.41 


17.24 


136.9 ± 7.4 


0.85 




TWIN 


144.7 ±6.0 


2.29 ±0.27 




132.5 ±5.7 


1.14 



Notes. — Columns are the same as in Table p] except for Col. (2); here, we list the instrument used for the observation. 



FRs), in which the rotation velocity is the depen- 
dent variable (y = logl/rot) and the photometric 
quantity is the independent variable {x = O = 

{MA,Msyn,logM,,Boll,logM.,MPA,logAfbar}). FitS tO 

the different ITFRs will be presented in i|8] 

For each photometric quantity O, we model the 
fTFR as a linear relation with an intrinsic Gaussian 
scatter of width a in log Kot : 

log Kot = a + bx (Mx - Ma,p), (25) 
log Kot = a + bx {M,yn - Msyn.p). (26) 
log Kot = a+ (-2.56) X (logM, -logM,,p), (27) 

Here, a is the zero-point, b is the slope, and Op is the 
pivot value of the ITFR. Both a and a have units of 
log(kms~^) and b has units of log(kms~^) mag~^. We 
will refer to the slope of the ITFR as steeper (shallower) 
if b is more (less) negative. In most cases, we set the 
pivot value Op to be the weighted mean of the fit sam- 
ple. Doing so makes the covariance between the error in 
a and the error in b negligible (in other words, changing 
the value of b by its l-cr uncertainty does not change the 
best-fit value of a). 

For completeness, we also consider forward TFRs, 
in which the independent variable y = O and the de- 
pendent variable x = log Kot- Fits to the forward TFRs 
will be presented in §9.3[ 

For each photometric quantity, we model the for- 
ward TFR as a linear relation with an intrinsic Gaussian 
scatter of width (Tfwd in O: 

Mx = flfwd + 6fwd X log(Kot/Kot,p), (28) 

Msyn = flfwd + 6fwd X log(Kot/Kot,p)- (29) 

logM. = flfwd + (-0.4fefwd) X log(Kot/Kot,p), (30) 

Here, both the zero-point af„d and (Tfwd have 
units of mag, and the slope 6f„d has units of 
mag (log(kms"^))"\ To directly compare the results 
of the forward and inverse fits, we convert the forward 
fit parameters to the equivalent inverse fit parameters: 
flconv = logKot,p + (C'p — afwd)/6fwd, feconv = l/6fwd, and 
o"conv = 5"fwd/fefwd. Thcsc Converted parameters have the 
same units as a, b, and a, as defined in Eqs. |25H27| 



Note that even in the ideal case in which there is 
zero uncertainty in the data, fits to the forward and 
inverse TFRs will not necessarily yield the same rela- 
tion. Moreover, fits to the forward TFRs are subject to 
Malmquist bias due to our sample selection cuts on ab- 
solute magnitude (in contrast, fits to the ITFRs are not 
affected by this bias). The inverse and forward fit direc- 
tions assume different underlying models and must be 
interpreted accordingly. In most classical TFR studies, 
forward fits were used to obtain distance measurements, 
based on the luminosities derived from the disk rotation 
velocities. In this work, we use fits to the ITFRs to ob- 
tain photometric estimators of disk rotation velocities. 

7.2 Fitting the TFRs 

We perform weighted maximum likelihood fits to the 
inverse and forward TFRs, defined by Eqs. |25f[27| and 
|28H30[ respectively. First, we describe how we assign 
weights to each galaxy. Then, we define the likelihood 
function used to determine the fit parameters and the 
bootstrap resampling method used to estimate their un- 
certainties. Finally, we describe the two-step method for 
fitting the TFRs for synthetic magnitudes Msyn. 

We assign to each galaxy a weight according to its 
stellar mass, so that the effective stellar mass function 
(SMF) of the fit sample matches that of the parent disk 
sample. Weights are given by the nearest integer to 100 x 
Wsamp(M*,Bcii), where 

Apparent (M,) iVsamp 



ap(M.) 



-^parent -^samp (-^'^* ) 



(31) 



Here, A'^parcnt ~ 169 563 galaxies, A'^samp is the total 
number of galaxies in the fit sample, and Npa,icnt{M,) 
and A'^samp {M, ) are the number of galaxies in the stellar 
mass bin containing M* , for the parent disk sample and 
the fit sample, respectively. We use logarithmic stellar 
mass bins of width « 0.21 dex. 

We find that the child disk sample SMF deviates 
from the parent disk sample SMF only at the lowest and 
highest stellar mass bins, with uigamp ~ 2 for galaxies 



with stellar masses below lO^ '^Mp) and 



1/2 
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Figure 19. Comparison between our measurements and 
those from P07 for 157 galaxies in the P07 sample (filled 
circles). Panels show A(logX), where X is as labelled, from 
top to bottom: rotation curve fit parameters -RtOi K: obsi 
inclination corrected rotation velocity Vgo, and inclination 
corrections sin 6. The error bars shown are our derived l-cr 
measurement uncertainties. Galaxies for which lAVcobsl is 
greater than 2a-{Vc obs) E'i'e marked by red circles in all pan- 
els (19 galaxies, or 12 per cent of the sample). Horizontal red 
solid lines mark mean values, (A(logX)); black dashed lines 
mark the zero axis. 



where K is a constant. 

We determine the three best-fit parameters, a, b, 
and a, by maximizing Eq. |32] using the Levenberg- 
Marquardt routine mpfit2dfuii in IDL (Markwardt 
20091. Then, we calculate the 1-a uncertainty in each 



parameter using bootstrap resampling. We generate 500 
sample realizations by randomly drawing A'' galaxies 
with replacement. We perform fits to each realization, 
and check that the distributions in the bootstrap pa- 
rameters are approximately Gaussian. 

To fit for the TFRs for synthetic magnitudes Mgyn, 
we first have to fix the coefficient a in its definition (c.f. 
Eq. |8|. To do this, we adopt a four-parameter model 

log Kot = 02 + 62 [Mx - Mx,p2 +ctx (colour)], (33) 

with a Gaussian intrinsic scatter of width a2 in log Kot- 
The likelihood function for this model is defined anal- 
ogously to Eq. |32[ Minimizing the likelihood function 
yields the best-fit coefficient a, and bootstrap resam- 
pling yields its 1-a uncertainty. Finally, the fits to the 
inverse and forward Msyn TFRs are determined from 
Eqs. |26| and |29| respectively, in a similar manner as for 
Mx, as described above. 

We present resuhs of weighted TFR fits in !|8] and 
i|9] In practice, though, we find that unweighted fits 
yield almost identical results. This is not surprising be- 
cause, as noted above, the SMFs of the child and parent 
disk samples are not significantly different. Therefore, 
for simplicity, we have not used weights in fits for the 
coefficient a in Mgyn, nor in the calculations performed 
in subsequent sections §10|and §11| 



7.3 Interpreting the TFR scatter 

Our TFR models include a Gaussian intrinsic scatter of 
width cr. In addition, we define the measured scatter, 
o"meas, to be the RMS (root-mean-square) of the TFR 
residuals {Ay)i = j/i — [a + fe x (xi — Xp)]. In the simplest 



for stellar masses above In contrast, the P07 

disk sample SMF is significantly different from that of 
the parent disk sample, especially for the highest stellar 



mass bins where Wb: 



0.1. 



The full dataset, with repeats, consists of a total of 
N = round(lOOwj) pairs of measurements, Xi 

and iji, with measurement errors ax,i and Uy^i. Assum- 
ing that these errors are uncorrelated, and that the mea- 
surements are normally distributed around their true 
values, X and y, with Gaussian widths given by their 
respective measurement errors, then the log-likelihood 
function is given by 



ln£ 



1 ~ 



2 , 2 



I 1,2 2 



equal for all galaxies, and in which our assumptions re- 
garding the Gaussian distributions for Xi and iji hold, 
o"mcas ~ 5-^ + Ty^i + {bax,i)^- In other words, the mea- 
sured scatter can be attributed partly to the measure- 
ment uncertainties for individual galaxies, and partly 
to the intrinsic scatter in the ITFR itself. Note, though, 
that what we have defined to be the "intrinsic scatter" in 
the TFR, cr, actually includes contributions from other 
sources of observational scatter that we have not for- 
mally accounted for in the measurement uncertainties, 
including errors in inclinations, slit PA misalignments, 
and intrinsic disk ellipticities. 

As noted in §6.5. 1| ellipticities in the potential in 
the plane of the disk e$ lead to noncircular motions 
and differences between the kinematic and photomet- 
ric major axis PAs. They also lead to a dependence of 
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the inclination correction on the viewing ang ler Ob- 
servations indicate that e$ « 0.1 ( Ryden|[2006 [Padilla 



fc Strauss||2008 and references therein), which in turn 
translate into an expected scatter in the (forward) TFR 
of 0.46 mag, even if kinematic major axis PAs were used 
( Franx fc de Zeeuw|1992 I . The expected scatter will be 
larger if photometric major axis PAs were used, as is 
the case in most TFR studies, including ours. In §12. 2| 
we compare the observed scatter in the TFR with the 
expected contributions from these systematic effects. 



8 CALIBRATED ITFRS 

In this section, we present weighted fits to the ITFRs 
for the 189 galaxies in the child disk sample using differ- 
ent photometric quantities O = {M\, M^yn, logAf*,Bcii, 
logMt^MPA, and logMbar} (all defined inj5| and rota- 
tion velocity amplitudes Vgo (defined in \ 6.4 1. These cal- 



ibrated ITFRs provide photometric estimators of disk 
rotation velocities applicable to the parent disk sample. 
The procedure for obtaining these fits were described 
in the previous section, !|7] Alternative fits will be pre- 
sented in the next section, iJ9l 



8.1 Ma ITFRs 

Table [8] lists the best-fit parameters a, b, and a for 
the calibrated Mx ITFR, for five SDSS bands, from 
weighted fits to the child disk sample (189 galaxies). 
Also listed are the pivot value xq = Mx^p, which has 
been set equal to the weighted mean of Mx, and the 
measured scatter CTmeas (defined in S7.3I. 

We confirm previous studies and find that the 



amount of scatter in the ITFR (both intrinsic and mea- 
sured) systematically decreases toward redder bands. 
Among the SDSS bands, we choose the i band to be 
the optimal one, as it yields a tighter ITFR relation 
than the bluer bands (in addition to being less affected 
by dust extinction). The scatter in the z band ITFR is 
similar, but this band has noisier photometry. 

8.2 Msyn ITFRs 

Table [9] lists the best-fit parameters for the calibrated 
Msyn ITFRs. We find that regardless of the choice of 
absolute magnitude band, the amount of intrinsic and 
measured scatter for the M^yn ITFRs is the same (0.035 
and 0.056 dex in logVgo, respectively). The zero-points 
and slopes are consistent with one another as well. We 
find that the Msyn ITFRs have smaller intrinsic scatter 
than the AIx ITFRS, by at least 2a; their measured scat- 
ter are also smaller. We conclude that synthetic magni- 
tudes are better photometric estimators of disk rotation 
velocities than single-band optical absolute magnitudes. 

If e<j. "C 1, the true inclination correction is sin 9(1 — 
e^cos2(/)) instead of sin 6, for a viewing ang le {0,<f>) l |Franx| 



Table[9]also shows the best-fit coefficients a (col. 4) 
in the definition of synthetic magnitudes Mgyn = Mx + 
a X (g ~ r). These are to be compared with the coef- 
ficients based on the colour-based estimates of Mt/Lx 
from BeU et al. (2003), q^'=" (col. 3). For each band A, we 
find that the best-fit a is within la of a^°". These colour 
corrections extrapolate toward a redder band, so that 
the synthetic magnitudes approximate stellar masses. 
We find similar results when using colours other than 
g — r. These results suggest that stellar mass is indeed 
a more fundamental parameter than luminosity in the 
TFR. Moreover, they provide independent support for 
the reliability of these colour-based stellar mass-to-light 
ratios (Kannappan et al. 



2002 



came to similar conclu- 
sions regarding the earlier Bell fc de Jong 2001 stellar 
mass-to-light ratios). 



8.3 M, 



,Bell) 



M, 



I, and Mbar ITFRs 



Table [To] lists the best-fit parameters for calibrated IT- 
FRs for stellar mass estimates, M,3eii and M«,mpa, and 
for baryonic mass estimates, Mbar- As expected from the 
systematic difference between M*,Beii and Mt,MPA (see 
Fig. |6] in S 5.4.3 1, the former has a pivot value higher 
by 0.15 dex, corresponding to a lower normalization in 
logKot, of 0.04 dex (at fixed M.). The M.,bo11 ITFR 
is steeper than the A'U^mpa ITFR at > Aa level. Their 
intrinsic and measured scatter are similar to each other 
and to the Afsvn ITFR. These results confirm our ear- 



lier conclusion (in [S.2) that colour-based stellar mass 



estimates are a good tracer of disk rotation velocity. 

Compared to the M.^eeii and M«_mpa ITFRs, the 
Mbar ITFR has larger intrinsic and measured scatter. 
True baryonic mass is expected to trace the disk ro- 
tation velocity more faithfully than stellar mass alone, 
since the both stellar and gas mass contributes to the to- 
tal disk potential. However, we do not find this to be the 
case for these baryonic mass estimates. Moreover, based 
on the increasing trend in gas mass fraction as stellar 
mass decreases, one would expect the baryonic ITFR to 
have a shallower slope than the stellar mass ITFR; con- 
trary to this expectation, we find a best-fit slope that 
is steeper by ~3cr. We conclude that Mbar is not an 
optimal photometric estimator for disk rotation veloci- 
ties. In other words, the relation between u — r colour 
and gas-to-stellar mass ratios (Eq. 141, on which our 



: de Zeeuwl 



|T992 i 



baryonic mass estimates are based, is not tight enough 
for our purposes (unlike the colour- Mt/L relations on 
which our stellar mass estimates are based). 



8.4 Comparison of Mi, M^yn, and A/,,bc11 ITFRs 

In this subsection, we compare three good photomet- 
ric estimators of disk rotation velocity, choosing one of 
each type we consider, namely: Mi, Msyn (i, 3 — r), and 
M«_Bcii. Table [TT] puts together the fit results from Ta- 
bles [8| jl0[ for easier comparisons. Figure [20] shows the 
observed relations (filled circles with 1-a error bars), the 
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Table 8. Calibrated Mx ITFRs for the child disk sample (189 galaxies) with Vrot = Vso- 



y 


X 


xo 


a 


b 


a 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


log Vso 


Mu 


-18.731 


2.188(0.026) 


-0.080(0.009) 


0.114(0.011) 


0.114 




Mg 


-19.903 


2.142(0.006) 


-0.129(0.008) 


0.068 (0.008) 


0.082 




Mr 


-20.375 


2.142(0.005) 


-0.130(0.007) 


0.056 (0.007) 


0.071 




Mi 


-20.558 


2.142(0.005) 


-0.128(0.006) 


0.049 (0.007) 


0.066 




M^ 


-20.649 


2.142(0.005) 


-0.119(0.005) 


0.049 (0.007) 


0.065 



Notes. — Col. (1): y is the dependent variable of the ITFR, or log Vrot- Col. (2): x is the independent variable of the ITFR, 
or C Col. (3): Xp is the pivot value for x, equ al to th e weighted mean of the fit sample. Cols. (4-6): Best-fit ITFR parameters 
and their l-tr uncertainties, as defined in Eqs. |25||26| a is the zero-point, in units of log(kms~^), b is the slope in units of 
log(kms~"'^) mag^-"^, and a is the intrinsic Gaussian scatter in units of log(kms~'^). Col. (7): CTmcas is the measured scatter in 
the ITFR, defined to be the RMS of the ITFR residuals (Aj/), = yi - [a + b X {x - Xp)]. 



Table 9. Calibrated Msyn ITFRs for the child disk sample (189 galaxies) with Viot = Vgo, where Mgyn = + a X {g — r). 



A 


colour 




a 




a 


b 


(J 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


9 


g-r 


-3.80 


-4.09 (0.91) 


-22.172 


2.159(0.004) 


-0.105(0.004) 


0.034(0.006) 


0.056 


r 


g-r 


-2.74 


-3.25 (0.72) 


-22.211 


2.160(0.004) 


-0.105(0.004) 


0.035 (0.005) 


0.056 


i 


g-r 


-2.16 


-2.42 (0.86) 


-21.965 


2.160(0.004) 


-0.107(0.004) 


0.035 (0.006) 


0.056 


z 


g-r 


-1.72 


-2.28(1.04) 


-21.987 


2.160 (0.004) 


-0.101 (0.004) 


0.036 (0.006) 


0.057 



Notes. — Col. (l)-(2): Band A and colour in the synthetic magnitude Mgyn = Mx + a X (colour). Col. (3): o^"" is the 
coefficient of the colour correction in the Bell et al. (2003) stellar mass-to-light ratio estimates for this band and colour. Col. 
(4): Best-fit coefficient a and its 1-cr uncertainty. Col. (5): Xp is the pivot value for x = M^yn. Col. (6-8): Best-fit ITFR 
parameters and their 1-cr uncertainties: a is the zero-point, in units of log(kms~^), b is the slope in units of 
log(kms~^) mag~^, and a is the intrinsic Gaussian scatter in units of log(kms~^). Col. (9): cr^neas is the measured scatter in 
the ITFR, defined to be the RMS of the ITFR residuals (Ay)^ = yi - [a + b X {x - Xp)]. 



best-fit ITFRs (thick solid lines), and tiie amounts of in- 
trinsic scatter in tlie relations (upper and lower dashed 
lines are displaced by ±la from the mean relation). Note 
that Msyn(i,ff - r) = Mi- 2.42(g - r) (c.f. Tablejof. 

We find that adding colour information to luminos- 
ity, in Msyn and Mt^cii, reduces the amount of scat- 
ter in the ITFR. In Figure [21] we directly compare the 
Mi and Msyn ITFRs (black inverted triangles and lines 
vs. red upright triangles and lines, respectively), with 
each pair of data points connected galaxy-by-galaxy. 
The flattening of the slope from one relation to the 
other is expected from the color-luminosity relation: 
brighter galaxies are redder and therefore have larger 
color corrections than fainter galaxies. The reduction in 
the scatter can also be explained similarly: outliers that 
lie well above the Mi ITFR relation (in particular, those 
with log Vso ^ 2.3) tend to be redder and have larger 
color corrections, putting them closer to the mean Msyn 
ITFR. 

Among these photometric quantities, we choose 
M,.Bcii to be the optimal photometric estimator of disk 
rotation velocity. It yields the one of the tightest rela- 
tions among those we have considered. It has a natural 
physical interpretation, as an estimate of stellar mass. 
Finally, it is straightforward to determine from photom- 
etry (as a simple linear combination of absolute magni- 
tude and color). 

Before we end this section, we study the distribu- 



tion of residuals from these ITFRs. Visual inspection 
of Fig. [20] indicates that the amount of scatter in the 
ITFRs varies over the full range in O. In particular, 
the relation seems to be tightest for intermediate values 
of O. To investigate further, we study the normalized 
residuals A(log Kot)i/o"T,i, where the numerator is the 
offset of a galaxy's measured rotation velocity from the 
best-fit ITFR, and the denominator is the "total scat- 
ter" associated with that galaxy, defined as a combina- 
tion of its measurement errors and the intrinsic scatter 
for the full ensemble, 

o"T,i = + cr^j + {ba^,if. (34) 

We find that the contribution from the uncertainty in 
the rotation velocity is the dominant term: for the child 
disk sample, (cr^)^''^ = 0.063 and (a^y^'^ = 0.081, 
0.073, and 0.073, for the M,, Msyn(i,5-r), and M,,bo11 
ITFRs, respectively. The contribution from ax,i is al- 
ways negligible because it is heavily downweighted by 
the coefficient « 0.01. 

Figure [22] shows distributions of the normalized 
residuals from the Mi, AIsyn{i, g — r), and M,,bg11 ITFRs 
(histograms in the top, middle, and bottom panels, re- 
spectively). For comparison, Gaussian distributions of 
unit width and of width equal to the standard devia- 
tion (s.d.) of the distribution are also shown (dashed and 
solid curves, respectively). The leftmost panels show dis- 
tributions for the full child disk sample and the other 
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Figure 20. ITFRs between rotation velocities Vgo and photometric quantities 0={Mi, log Af,, and Meynii, g — r)} (left, middle, 
and right panels, respectively). Observations for the 189 galaxies in the child disk sample are shown as filled circles with 1-a 
error bars. Thick solid lines show the best-fit ITFRs and dashed lines indicate the amount of intrinsic scatter (they are displaced 
by ±5" from the mean relations). Vertical dotted lines divide the sample into three bins of roughly equal number. Labels list 
the pivot values Op, best-fit parameters — zero-point a, slope b, intrinsic scatter a — of the ITFR, and measured scatter crmoas 
in log Vgo . 



Table 10. Calibrated A/, Bclh Af,.MPAi and M^ar ITFRs for the child disk sample (189 galaxies) with Vrot = Vso- 



y 


X 




a 


b 


u 


C mcas 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


log Vgo 




10.102 


2.142(0.004) 


-0.111 (0.004) 


0.036 (0.005) 


0.056 




A^«,MPA 


9.952 


2.140(0.004) 


-0.093 (0.004) 


0.034 (0.006) 


0.056 






10.355 


2.144(0.005) 


-0.127(0.006) 


0.053 (0.007) 


0.071 



Notes. — Columns are the same as in Table |8] 



panels show distributions for the faint (low-Af*), inter- 
mediate, and bright (high-M«) bins (from left to right). 

We expect that the distributions for the full sample 
would be Gaussians of unit width. We find that the dis- 
tribution widths are indeed close to unity, although the 
distributions are more centrally-peaked than a Gaus- 
sian (leftmost panels in Fig. 221. We find that the dis- 
tribution widths vary with stellar mass and luminosity. 
The intermediate bins have the narrowest distributions, 
with s.d. values less than unity. On the other hand, the 
bright /high- Af, bins have the widest distributions, with 
s.d. ~ 1.3. These galaxies have better S/N measure- 
ments, and therefore tend to have smaller measurement 
errors on average (with {<JtY^^ = 0.05, compared to the 
sample mean of 0.08). The faint/low-Af* bins have the 
largest absolute residuals, but those galaxies also have 
the largest measurement errors. 

It is interesting to compare the distributions for the 
Mi ITFR on the one hand, and the Af. and A//syn IT- 



FRs on the other. We find that the reduction in scatter 
when going from Mi to Af* (or Msy-n) can be mostly 
attributed to galaxies in the middle bin, with the s.d. of 
the distribution in normalized residuals dropping from 
0.8 for Mi to 0.7 for Af^yn (or Af,,Beii). 



9 ALTERNATIVE TFR FITS 

9.1 Alternative Kot amplitudes 

In i|8] we have presented ITFR fits with Kot ~ Vso, 
our default rotation velocity amplitude definition. Now, 
we present ITFR fits with alternative rotation velocities 
V2.2 (the rotation velocity evaluated at a radius 2.2Ra) 
and Vc (the asymptotic rotation velocity of the best-fit 
arctangent model). Table 



12 



shows the best-fit param- 
eters for the Mi, Aft,Bcii, and Msyn{i,g — r) ITFRs for 
Viot = V2.2 and Vc- We find that a is robust to the choice 
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Table 11. Calibrated Mi, Msyn{i, g — r), and Af, Boll ITFRs for the child disk sample (189 galaxies) with V^ot = Vgo (compiled 
from Tables [spoil. 



y 


X 


Xp 


a 


b 


a 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


log Vgo 




-20.558 


2.142(0.005) 


-0.128(0.006) 


0.049 (0.007) 


0.066 




Msyn{i,g - r) 


-21.965 


2.160(0.004) 


-0.107(0.004) 


0.035 (0.006) 


0.056 




'M"*,Bcll 


10.102 


2.142(0.004) 


-0.111 (0.004) 


0.036 (0.005) 


0.056 



Notes. — Columns are the same as in Table [s] 



V 




A 





2.5 
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„2.1 
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Figure 21. Comparison between Mi and Mgyn ITFRs (black 
inverted triangles and lines vs. red upright triangles and lines, 
respectively). Pairs of points for every galaxy are connected 
by gray horizontal lines. The amount of intrinsic scatter in 
the ITFRs is indicated by the dashed lines, which are dis- 
placed by it 5- from the mean relations. Vertical dotted lines 




divide the sample into three equal bins in Mi and 
(black and red dotted lines, respectively). 



Ms, 



of rotation velocity definition, and consistent with the 
Bell et al. (2003) coefficient q''°" = -2.16 (as is the case 
for Vrot = Vgo). 

The panels of Figure [23| compare the M*,bc11 ITFR 
with Vrot = Vso (red solid and dashed lines, on both 
left and right panels), with the observed ITFRs with 
Kot = V2.2 and Vc (black solid and dashed lines in the 
left and right panels, respectively). As expected from the 
systematic difference between the V2.2 and Vso, the zero- 
point of the V2.2 ITFR is 0.02 dex (5 per cent) lower than 
that of the Vso ITFR. The slope of the ITFR is slightly 
steeper for the former, but only by ~ la. The amounts 
of intrinsic and measured scatter are similar for the two 
relations. These results confirm our earlier assertion (in 
S6.4I that the ITFR fits are relatively insensitive to the 



choice of rotation velocity amplitude definition, for rea- 
sonable choices for the optical (i.e., evaluation) radius. 
In contrast to Vso and V2.2, using K as the rotation 



velocity involves extrapolation well beyond the observed 
rotation curve. This can lead to systematic overestima- 
tion of the disk rotation velocity, especially in galax- 
ies with large turn-over radii (i.e., slowly-rising rotation 
curves). We find that the best-fit Vc ITFR has a sig- 
nificantly higher zero-point, shallower slope, and larger 
intrinsic scatter than the Vgo ITFR. The right panel 
of Fig. [23] shows that the two relations coincide at the 
high-mass end. This is consistent with our observation 
that the highest M« galaxies have small turn-over radii 
(i.e., steeply-rising rotation curves), and consequently 
Vc/Vso ratios close to unity (see Fig. 12 1. 



9.2 Alternative internal extinction corrections 

We also present fits to the A'I\ ITFRs for alternative 
internal extinction corrections. Our default internal ex- 
tinction corrections are applied to the mean orientation 
of the child disk sample, and chosen so as to minimize 
the s ize, and th eref ore uncertainty, of the corrections 
(see ijirl. Table 



13 



hsts ITFR fit parameters for Ml~", 
with extinction corrections applied to face-on orienta- 
tion (usually adopted in the literature, in particular by 
P07), as well as , with no extinction corrections at 
aU. 

For easy comparisons, we have chosen the pivot val- 
ues Ma,p to be that used for the default fit. As expected 
from the differences in the amount of extinction correc- 
tions, the normalization of the ITFR is lower by ~0.04 
dex for M^-°, and higher by ~ 0.003 for M^^. The 
amount of measured scatter is similar in all three cases; 
the amount of intrinsic scatter is slightly smaller for 
the extinction-corrected magnitudes than for the uncor- 
rected ones (because of the additional uncertainty from 
the extinction corrections). 

As a tes t of o ur internal extinction corrections, we 
will show in S 10.1 that the velocity residuals in the Mf"'" 



ITFR are correlated with disk axis ratios, while those 
for the Mi ITFR are not. 



9.3 Alternative fit directions 

We present calibrated forward fits to the Mi, Msyn{i, g — 

In these fits, the pho- 



14 



r), and M,,bc11 TFRs in Table 
tometric quantity O serves as the dependent variable, 
and Kot = Vso serves as the independent variable (see 
Eqs. [2ip0l ). We list the best-fit forward TFR param- 
eters Ofwd, fcfwd, and 5"fwd, as well as the converted fit 
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all galaxies faint/low-mass bin middle bin bright/high-mass bin 



0.5 

'0.4 

0.3 

0.2 

0.1 

0.0 
0.5 

■ 0.4 

0.3 

0.2 

0.1 

0.0 
0.5 

■ 0.4 
0.3 
0.2 
0.1 
0.0 



Sid. = 0.990 
= 0.081 



Mi 




Sid. = 0.835 
(4);^^ = 0.108 



Sid. = 0.808 
(4)f/2 = 0.067 



Sid. = 1.253 
((T|,)f/2 = 0.059 




o 



o 

"o 
c 
o 

o 
CO 



sid. = 0.987 
(4)f/2 = 0.073 



Sid. = 0.901 - 

(4)1/2 ^ g Qgg 



Sid. = 0.709 - 
(4)f''2 = 0.082 



d. = 1.253 
/2 = 0.051 



M 



syn 




Sid. = 0.989 

(4)1/2 ^ Q 



Sid. = 0.904 - 
(4)f/2 = 0.084 



Sid. = 0.677 - 
(4)1/2 = 0.081 







Sid. = 1.286 - 
(4)1/2 = 0.050 



-3-2-10 1 2 3 
A(log V8o)/CTr 
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Figure 22. Distribution of normalized residuals A(log V8o)/o't for the Mi, Msyn(i, <? — r), and A/, ITFRs (top, middle, 
and bottom panels, respectively). For each galaxy, a^^i = (S^ + o-'y ^ -\- ^^ct^ j)l/^ is the total scatter, defined to be the root- 
square-sum of the intrinsic scatter in the ITFR and the galaxy's measurement error. We show distributions for the full child 
disk sample (left-most panels), as well as for the sample split into three bins in Mi, Msyn, and M,, Bell (rest of the panels; from 
left to right). Each panel is labeled by the standard deviation of the distribution (s.d.) and the RMS value of cry. A Gaussian of 
width equal to the measured s.d. is shown by a solid curve, and one with unit width is shown by a dashed curve, for comparison. 



Table 12. Calibrated Mi 
and Vc- 



Afsyn(*,9 — r), and Af, bcII ITFR, fits for the child disk sample (189 galaxies) with Vj-ot = Vgo, V2.2, 





y 


X 




a 


b 


(J 


fj 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


lo 


5^80 


Mi 


-20.558 


2.142(0.005) 


-0.128 (0.006) 


0.049 (0.007) 


0.066 






Msyn{i,g-r) 


-21.965 


2.160(0.004) 


-0.107 (0.004) 


0.035(0.006) 


0.056 






M.,Bell 


10.102 


2.142 (0.004) 


-0.111 (0.004) 


0.036 (0.005) 


0.056 


lo 




M, 


-20.558 


2.120 (0.004) 


-0.133 (0.006) 


0.052(0.007) 


0.065 






Msy-n{i,g-r) 


-21.799 


2.142(0.004) 


-0.107(0.004) 


0.035(0.005) 


0.056 






Af*,Bcll 


10.102 


2.120(0.004) 


-0.116(0.004) 


0.036 (0.005) 


0.055 


lo 




M, 


-20.558 


2.205(0.005) 


-0.112(0.006) 


0.053 (0.007) 


0.093 






Msyn(i, g - r) 


-21.799 


2.142 (0.004) 


-0.107 (0.004) 


0.035(0.005) 


0.056 






M..B0II 


10.102 


2.205 (0.004) 


-0.095 (0.004) 


0.046 (0.006) 


0.087 



Notes. — Columns are the same as in Table |8] Results for the y = logVso are the same as in Table [TT] 



parameters a, 



conv , i-'conv 



in S7.ll 



and CTconv (defined after Eq. 30 



Compared with the inverse fits (shown in the lower 
half of Table 14 for easy comparisons), the forward fits 
have slightly larger intrinsic scatter (by ~lcr), and sig- 
nificantly steeper slopes (feconv < h). Figure 24 shows 



the inverse and forward fits to the A/*3cii TFR (black 
solid and dotted lines, respectively). 

The steeper slopes from the forward fits can be (at 
least, partly) attributed to a Malmquist-type bias due to 
our applied absolute magnitude cuts. Based on a Monte- 
Carlo experiment performed by P07 (see their Appendix 
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log M..Bell (Mo) log M,,Bell (Mq) 



Figure 23. Calibrated Bell ITFRs with rotation velocities V2.2 and Vc (left and right panels, respectively) for the child disk 
sample (filled circles with l-cr error bars). In each panels, the best-fit ITER is shown by the black solid line and the amount 
of intrinsic scatter is shown by the dashed lines (displaced by ±5- from the mean relation). For comparison, the Af* bcII ITFR 
with Vrot = Vgo is shown by red solid and dashed lines. 



Table 13. Calibrated Mx ITFR fits for alternative internal extinction corrections for the child disk sample (189 galaxies). 



y 

(1) 


X Xp a b 
(2) (3) (4) (5) 


(6) 






C mcas 

(7) 




With default internal extinction corrections (same 


as Table 


8 





logVso Mu -18.731 2.188(0.026) -0.080(0.009) 0.114(0.011) 

— Mg -19.903 2.142(0.006) -0.129(0.008) 0.068(0.008) 
-— Mr -20.375 2.142(0.005) -0.130(0.007) 0.056(0.007) 

— Mi -20.558 2.142(0.005) -0.128(0.006) 0.049(0.007) 
— M^ -20.649 2.142(0.005) -0.119(0.005) 0.049(0.007) 



With internal extinction corrections to face-on orientation 



log Vso 


Mr" 


-18.731 


2.142 (0.034) 


-0.072(0.015) 


0.112(0.016) 


0.112 




Mr" 


-19.903 


2.079 (0.007) 


-0.117(0.007) 


0.064(0.009) 


0.082 




Mr" 


-20.375 


2.087(0.006) 


-0.120(0.006) 


0.052(0.008) 


0.071 




Mr" 


-20.558 


2.098 (0.005) 


-0.119(0.005) 


0.046 (0.008) 


0.066 




Mr" 


-20.649 


2.111 (0.005) 


-0.112(0.004) 


0.047(0.007) 


0.065 


Without internal extinction corrections 


log V8O 


M^<- 


-18.731 


2.190 (0.022) 


-0.081 (0.011) 


0.114(0.017) 


0.114 




MNC 


-19.903 


2.147 (0.006) 


-0.123(0.008) 


0.072(0.009) 


0.083 




MNC 


-20.375 


2.147 (0.005) 


-0.126(0.006) 


0.060(0.009) 


0.073 




MNC 


-20.558 


2.145 (0.005) 


-0.127(0.005) 


0.052(0.008) 


0.066 




MNC 


-20.649 


2.145 (0.005) 


-0.119(0.005) 


0.049 (0.008) 


0.065 



Notes. — Columns are the same as in Tabled 



B), the effect of the cuts is to steepen the slope and re- 
duce the intrinsic scatter. We do not attempt to correct 
for this bias and present the forward fit results at face 
value. Throughout the rest of the paper, we focus on 
inverse fits, which are not affected by this bias. 

We have also performed orthogonal fits to the 
TFRs, in which offsets are measured in the direction 
orthogonal to the best-fit relation itself (rather than 
in the horizontal or vertical directions). We find that 



the best-fit parameters from the orthogonal fits are al- 
most identical to those for the inverse fits. This is not 
surprising because the total measurement uncertainty 
is dominated by the rotation velocity uncertainty, i.e., 
ai(logl/rot) is much larger than 6 x ai{0) (c.f. Eq. 
and discussion afterwards). 



34 
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Table 14. Calibrated forward fits to the TFRs for the child disk sample (189 galaxies) with independent variable x = logVgo 
and dependent variable y = O = {Mi, Msyn{i, g — r), log Af* bcIi}- 



o 

(1) 


Xp 

(2) 


Ifwd 

(3) 


(4) (5) 


^conv 

(6) 


0"fwd 
(7) 


f^conv 

(8) 


Forward fits to the TFR 


Mi 

Msyn 
M.,Bcll 


2.142 
2.160 
2.142 


-20.602 (0.033) 
-21.983(0.032) 
10.054(0.033) 


2.131(0.006) -5.807(0.307) - 
2.115(0.004) -7.716(0.310) - 
2.133(0.004) -7.465(0.281) - 


0.172 (0.009) 
0.130 (0.005) 
0.134 (0.005) 


0.344 (0.041) 
0.312(0.046) 
0.307(0.043) 


0.059 (0.007) 
0.041 (0.006) 
0.041 (0.006) 






Inverse fits to the TFR, with Viot = Vgo (same 


as in Table 11 


1 


M, 

Msyn 
M,,Bcll 


-20.558 
-21.965 
10.102 




2.142(0.005) ... - 
2.160(0.004) ... - 
2.142(0.004) ... - 


0.128 (0.006) 
0.107 (0.004) 
0.111 (0.004) 




0.049 (0.007) 
0.035(0.006) 
0.036 (0.005) 



Notes. — Col. (1): O is the photometric quantity — the dependent (independent) variable in the forward (inverse) fits. Col. 

(2) : Xp is the pivot value for the independent variable — (logV8o)p for the forward fits and Op for the inverse fits. Cols. 

(3) , (5), and (7): Best-fit forward TFR parameters — ifwdi ^fwdi and '^(wd — and their 1-cr uncertainties, as defined in 
Eqs. |28||30| Cols. (4), (6), and (8): Converted forward TFR parameters to be compared with the inverse fit parameters. 




1.0 9.5 10.0 10.5 11.0 

logM.,Bell(M0) 



11.5 



Figure 24. Comparison of alternative fits to the Bell 
TFR for the child disk sample (189 galaxies). Galaxies with 
D/T > 0.9 (137 galaxies) and D/T s£ 0.9 are shown by 
yellow and red circles, respectively (shown with l-cr error 
bars). Solid and dashed black lines show the best-fit ITFR 
and itlcr bounds; the dotted black line shows the best-fit 
forward TFR. Red dashed lines show ±15" bounds for the 
best-fit ITFR for galaxies with D/T > 0.9. 



9.4 Alternative sample: D/T > 0.9 

For various reasons, disk-dominated galaxies are ex- 
pected to form a tighter TFR than the general pop- 
ulation. The lack of a significant bulge component 
makes the measurement of various photometric quan- 
tities cleaner (i.e., less noisy). In addition, late-type spi- 
rals have been found to have smaller intrinsic disk el- 
lipticities than earlier types ( Ryden|2006 1, leading to a 
smaller contribution to the scatter due to random view- 
ing angles. Disk-dominated galaxies tend to be kinemat- 



ically homogeneous systems as well. Figure [24] shows 
that galaxies in the child disk sample with D/T > 0.9 
(135 galaxies; yellow circles) indeed form a tighter TFR 
than the full sample. 

We have performed weighted ITFR fits to the disk- 
dominated subsample of galaxies with D/T > 0.9. The 
red dashed lines in Fig. |24| mark ±15- from the best-fit 
relation. We find a = 0.027 ± 0.009 and ericas = 0.049, 
compared to a = 0.036 ± 0.005 and amcas ~ 0.056 for 



the full child disk sample (Table 111. We conclude that 



inclusion of only the disk-dominated galaxies yield a 
tighter ITFR, at the ~lo- level. 



9.5 Comparison with POT fits 

To make a fair comparison with the results presented 
in P07, we have performed ITFR fits to the P07 galaxy 
sample using the same photometric quantity and rota- 
tion velocity amplitude definition used by P07. The pho- 
tometric quantity used is absolute magnitude corrected 
for internal extinction to face-on orientation, M;~°, and 
the rotation velocity amplitude definition used is Vgo- 

First, we perform unweighted fits to P07 measure- 
ments (both for the x and y variables in the ITFR) 
and reproduce the results published in Table 4 of P07. 
Then, we perform weighted fits to our own measure- 
ments. Table [15] lists the three different sets of best- 
fit ITFR parameters. The relations derived using P07 
measurements and ours are consistent with one another 
within their reported intrinsic scatter, but the two de- 
viate slightly at the bright end (our measurements yield 
a best-fit ITFR with a shallower slope and lower nor- 
malization). The slight differences can be traced to dif- 
ferences in the determination of Vso (c.f. [ 6.5.3 1 and 
internal extinction corrections to Mi in the two meth- 
ods. 
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Table 15. Comparison of ITFR fits for the POT and child disk samples (with 157 and 189 galaxies, respectively). 





y X Xp a b 
(1) (2) (3) (4) (5) 


nicas 

(6) (7) 


ITFR parameters from Table 4 of P07 


lo 


gVso Ml'" -21.327 2.212(0.005) -0.130(0.005) 


0.061(0.005) 




Unweighted ITFR fits for the P07 galaxy sample, using 


P07 measurements 


lo 


gVso M^-" -21.327 2.207(0.005) -0.127(0.005) 


0.058(0.005) 0.067 




Calibrated ITFR fits for the P07 galaxy sample, using 


our measurements 


lo 


gVso M|"" -21.327 2.195(0.007) -0.122(0.006) 


0.055(0.009) 0.102 



Notes. — Columns are the same as in Table H 



10 TFR RESIDUAL CORRELATIONS 

In this section, we study correlations between velocity 
residuals from the calibrated ITFRs (derived in !|8] and 
summarized in Table [TT| and various galaxy properties, 
namely: disk axis ratios in §10. 1[ galaxy colours in §10. 2[ 
and disk sizes in §10. 3| 

We calculate two kinds of correlation coefficients: 
the Pearson linear correlation coefficient r and the 
Spearman rank correlation coefficient p. The former as- 
sesses how well the relationship between two variables 
can be described as linear {r — ±1 means there is a per- 
fect positive/negative linear correlation), while the lat- 
ter assesses how well the relationship can be described 
as monotonic, regardless of the form of the relationship. 
We also quote the two-tailed significance of the devia- 
tion of p from zero, Sig(p). A small value (i.e., smaller 
than 0.01) indicates a significant correlation. We also 
perform unweighted linear fits using a similar procedure 



as for the TFR fits (described in [ 7.2 1 and quote best-fit 
slopes b. 



10.1 Correlations with disk axis ratio 

To test our internal extinction corrections (defined in 
S 5.1 1, we study the correlation between disk axis ratios 
qd and velocity residuals from the Mf^ and Mi ITFRs, 



i.e., before and after internal extinction corrections (c.f. 
Table 13 1. We expect to see a positive correlation for 



the ITFR residuals, because more highly-inclined 

galaxies tend to be more heavily affected by dust and 
artificially displaced from the TFR toward fainter lumi- 
nosities. Hence, they tend to have larger (i.e., more posi- 
tive) velocity residuals. If our internal extinction correc- 
tions are valid on average, we expect the correlation to 
decrease or disappear altogether in the A4i ITFR resid- 
uals. 



Figure 25 shows log(l/gd) versus velocity residu- 
als A[log V8o(Mf<^)] and A [log Vso (Mi)] (left and right 
panels, respectively). The best-fit linear relations for 
179 (out of 189) galaxies in the child disk sample with 
qd < 0.6 are shown by solid lines on both panels. (Re- 
call that our target selection cut is based on giso, not qd; 
a number of child disk sample galaxies have qd slightly 
above 0.6.) For Mf"'", the best-fit slope is non-zero and 
positive, b = 0.11 ±0.06, while for Mi, the best-fit slope 



is consistent with zero, b = 0.02 ± 0.06. The correla- 
tion coefficients indicate a weak positive correlation in 
the former case, and a weak negative correlation in the 
latter, albeit neither correlation is very significant. We 
conclude that the applied internal extinction corrections 
are adequate on average. 



10.2 Correlations with galaxy colour 

The observed colours of a galaxy encapsulate proper- 
ties of its stellar population, in particular, its average 
stellar mass-to-light ratio M^/L. In !|8j we showed that 
by adding color information to the luminosity (i.e., us- 
ing M, or Msyn, instead of Mi), we can reduce the in- 
trinsic scatter in the ITFR by ~ 30 per cent (or 2a"). 
Now, we explicitly show that in doing so, we have suc- 
cessfully extracted most of the available information 
from galaxy colour to yield the tightest possible TFR. 
In other words, most of the correlation between galaxy 
colours and velocity residuals from the Mi ITFR is re- 
moved when Msyn is used instead. 

Figure[26]shows the correlation between g — r colour 
and A(log Vgo) for 189 galaxies in the child disk sample. 
Left and right panels show residuals from the best-fit Mi 
and Msyn ITFRs, respectively. Best-fit linear relations 
are shown as solid lines in both panels. Labels show 
the best-fit slope b and its 1-a uncertainty, the Pearson 
linear correlation coefficient r, the Spearman rank cor- 
relation p and the two-tailed significance of its deviation 
from zero. 

Clearly, we find a significant positive correlation be- 
tween galaxy colours and Mi ITFR residuals. The best- 
fit slope is 0.23 ± 0.04, and the Spearman rank correla- 
tion coefficient p = 0.39, with very high significance. On 
the other hand, we find no significant correlation for the 
Msyn ITFR residuals. The best-fit slope 0.04 ± 0.03 is 
consistent with zero, and p — 0.01. Indeed, the addition 
of colour information in A/syn removes essentially all of 
the residual correlation. 

Notably, the reddest galaxy in the sample 
SDSSJ203523-06 is an extreme outlier from the Mi 
TFR (with g-r = 1.31 ± 0.42 and A[log V8o(Mi)] = 
0.22), but lies within 25 of the M. and Msyn ITFRs 
(with log Vso residuals 0.03 and 0.05, respectively). This 
edge-on galaxy is an extreme case, with a prominent 
dust lane, disk axis ratio qd = 0.064 (smallest in the 
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Figure 25. Correlation between disk axis ratios and log Vgo residuals from the Mf-^ and Mi ITFRs (i.e., absolute magnitudes 
before and after internal extinction correction) (left and right panels, respectively). Filled circles are shown with l-cr error bars, 
for 179 (out of 189) galaxies in the child disk sample with < 0.6. Solid lines show the best-fit linear relations; horizontal 
dashed lines show the zero level. Labels show the best-fit slope b and its 1-a uncertainty, the Pearson linear correlation coefficient 
r, the Spearman rank correlation p and the two-tailed significance of its deviation from zero. 



sample) and M^/Li = 9.7 (largest in the sample; for 
comparison, the second largest value is 3.66, and the 
mean is 1.86), so it is gratifying to find that the color 
corrections seem to have worked in this case. 

Now, we show that the observed correlation be- 
tween g — r colours and Mi ITFR velocity residuals can 
be attributed to the variation in M,/L with colour. Red- 
der galaxies have larger Mt/L, so they tend to be fainter 
than average for a given rotation velocity. Assuming the 
BeU et. al. (2003) relation, M./L, oc (g - r)" ***^, we can 
predict the slope of the residual correlation to be 



A[logT/8o(MO] 
A(g - r) 



A [log 140 (Af.)] ^ A[log(Af./LO 



A(logM,) A{g-r) 
= (-2.5)(-0.111±0.004)(0.86) 
= 0.24 ±0.09, 



where we have used the best-fit slope of the M, ITFR 
from Table [To] This predicted slope is very close to the 
best -fit value (dotted vs. solid lines in the left panel of 
Fig.[26|If] 

Putting these results together, we form a fully con- 
sistent picture. The addition of colour information to 
Mi via M* and Msyn corrects for the variation in M^/ Li 
with colour. Doing so successfully removes most of the 
residual correlation with galaxy colour, and yields to 
significantly tighter ITFRs. 



10.3 Correlations with disk size 

Recent observational studies of the TFR have consis- 
tently found a lack of correlation between TFR residu- 



als and disk galaxy scale lengths (P07; [Courteau et al.| 



2007[ ). [Courteau fc Rix| ( [1999| ) noted that this lack of 
correlation argues that the disk mass does not con- 
tribute most of the disk rotation velocity at the opti- 
cal radius. However, Button et al. (2007) noted that 
there can be other interpretations as well. Meanwhile, 
[Gnedin et al.[ ( |2007[ ) emphasized that the observed lack 
of residual correlations plays a key role in their conclu- 
sions about the baryon fractions of disk galaxies. For 
a pure self-gravitating disk model, VJ-ot oc R^^ at fixed 
M«, so a strong negative correlation is expected with 
9 log Vrot /9 log 7?d ~ —0.5. On the other hand, a strong 
positive correlation is expected for a pure NFW dark 
matter halo model, with an expected slope of +0.5 in 
the inner regions. 

We study correlations between velocity residuals 
from the Af, ITFR and disk size offsets from the mean 
size-mass relation, A(log_Rd) = log(i?d/-Rd(A/,)). Fig- 
ure |27| shows the log A{f,-log iid relation for 189 galaxies 
in the child disk sample (filled circles), and the best-fit 
mean relation R^{M,) (solid line). Using a procedure 



similar to that used to fit the TFRs (c.f. i 7.2 1, we find 



If one instead assumes a pure gravitational disk model, in 
which V'j.Q^ (X Mi, at fixed scale length (as was done in P07), 
one would predict a different slope, equal to 0.86/2 = 0.43, 
which is inconsistent with the observations. 



log 



Rd 
kpc 



= (0.51 ±0.010) (0.20 ±0.02) X 
A/, 



log 



Mr, 



10.17 



(35) 
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Figure 26. Correlation between g — r colours and logVgo residuals from the Mi and Msyn ITFRs (left and right panels, 
respectively). Filled circles are shown with 1-a error bars for the 189 galaxies in the child disk sample. Solid lines show the 
best-fit linear relations; horizontal dashed lines show the zero level. The dotted line in the left panel has a slope of 0.24, the 
predicted slope from the Bell et al. (2003) scaling relation M^/L oc [g — r)"-^*. Labels show the best-fit slope b and its 1-a 
uncertainty, the Pearson linear correlation coefficient r, the Spearman rank correlation p and the two-tailed significance of its 
deviation from zero. 



with a best-fit Gaussian intrinsic scatter in logi?d of 
a = 0.14 ±0.01 dex. 

Figure [28] shows the correlation between velocity 
residuals A[log V8o(A^.)] and disk size ofTsets A(log J?d) 
for the 189 galaxies in the child disk sample. Confirming 
the results of prior studies, we find no evidence for a 
correlation: r = 0.05 and p — 0.04. The best-fit linear 
relation has a slope consistent with zero, b = 0.012 ± 
0.034 (solid line). 

Now, we go one step further and repeat the analysis 
for three bins in stellar mass (with 64, 64, and 61 galax- 
ies in the low, intermediate, and high stellar mass bins, 
respectively). Figure[29]shows the results for each stellar 
mass bin (as labeled; the legend here is similar to that 
in Fig. 28 1. Although the best-fit slopes are not close 



to either of the predictions from the pure disk and pure 
DM models (—0.5 and -1-0.5, respectively; dot-dashed 
lines), we find a decreasing trend in the best-fit slope 
(changing sign from slightly positive to slightly nega- 
tive) with increasing stellar mass: b = 0.087 ± 0.055, 
0.00 ±0.03, -0.04 ±0.08, for the low, intermediate, and 
high stellar mass bins, respectively. The correlations are 
weak, but reflect the same trend with increasing stellar 
mass: p = 0.25, —0.08, and —0.15, with corresponding 
Sig(p) = 0.05, 0.5, and 0.2, for the three bins, respec- 
tively. 

The observed trend indicates that the stellar mass 
(or baryon) fraction within the optical region of disk 
galaxies increases systematically with stellar mass, over 
the range of stellar masses we consider. In §11| we ex- 
plicitly calculate stellar mass fractions and confirm the 
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Figure 27. Relation between stellar masses bcII B.nd 
disk scale lengths for 189 galaxies in the child disk sam- 
ple (filled circles with l-cr error bars). The best-fit relation 
logiJ(j = a + 6(logMH, — 10.17) is shown by the solid line, 
and the dashed lines are displaced from the mean relation by 
ilCT, the best- fit Gaussian intrinsic scatter. Best-fit param- 
eters are listed in the upper left corner, together with their 
l-cr uncertainties. 
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Figure 28. Correlation between velocity residuals from the 
M« ITFR, A(logV8o), and disk size offsets A(log_Rd), de- 
fined relative to the mean relation log R^(Mi,) (given by 
Eq. \35\ . The best- fit linear relation has a slope consis- 
tent with zero (solid line). Predicted trends for a pure self- 
gravitating disk model (slope = —0.5) and a pure NFW DM 
halo model (slope = -1-0.5) are also shown (dot-dashed lines). 



trend with stellar mass suggested by these residual cor- 
relations. In §12. 4[ we discuss the interpretation of these 
results. 



11 DYNAMICAL-TO-STELLAR MASS 
RATIOS 

We calculate dynamical-to-stellar mass ratios within 
the optical radius Rso for the 189 galaxies in the 
child disk sample, denoted by (Mdyn/Mt)opt = 

(Mdyn/MH,)(i?80) = Mdyn(7?80)/M,(i?80) = 

Afdyn,opt/M*,opt. We adopt an empirical definition 
that depends straightforwardly on directly-observed 
quantities 



Mdyn 



VioRso/G 



opt 



0.8M,,Be 



[JC{D/T)s 



' kpc(kms 



(36) 



is the grav- 



where G = 4.3012 x 10' 
itational constant, IC = 1.34 is a geometrical factor that 
corrects for the flattened potential of the disk (assum- 
ing an exponential disk-|-NFW DM halo model; see Ap- 
pendix for the derivation), and {D/T)so is the disk-to- 
total mass ratio within _Rgo, 



iD/T)s 



1 



-flso/fld 



(1 + Rso/Rd) 



0.8 



x{D/T). (37) 



For simplicity, we have assumed that the average stel- 
lar mass-to-light ratio within _R80 is equal to its global 
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Figure 29. Correlations between velocity residuals from the 
M* ITFR, A(logV8o), and disk size offsets A(logi?d) for 
the child disk sample divided into three stellar mass bins, 
as labelled (with increasing M, from top to bottom). Solid 
lines show the best-fit linear relations, dashed horizontal lines 
show the zero level, and dot-dashed lines show the predicted 
trends for a pure self-gravitating disk model (slope = —0.5) 
and a pure NFW DM halo model (slope = -1-0.5). 



value, so that M*, opt/A/, « 0.8 by definition of Rso (ac- 
tually, the value of Mt,/L varies with radius, as indicated 
by color gradients along the disk). Note that VsqRso/G 
is the dynamical mass within _R80 for a spherical mass 
distribution with a circular velocity equal to I^80, and 
that we have assumed a fixed Kroupa IMF in the defi- 
nition of M*,Bell. 

We also calculate dynamical-to-baryonic mass ra- 
tios within Rso 
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dyn 



Mb 



dyn 



M, 
Mbai 



(38) 



where Mbar are baryonic mass estimates based on u — r 
color-based gas-to-stellar mass ratios from |Kannappanl 
(c.f. S5.5I. For simplicity, we have assumed that 



the average gas-to-stellar mass ratio within Rso is equal 
to its global value. Actually, the radial extent of HI is 



usually larger than the optical radius (Broeils & Rhee 



19971 



Figure [30| shows the relation between M«,opt and 
Mdyn.opt for the 189 galaxies in the child disk sample, 
using stellar mass estimates M.^cii and M,,mpa (top 
and bottom panels, respectively). The scatter in this 
relation is clearly being driven by disk size (red tri- 
angles, yellow circles, blue inverted triangles show the 
upper, middle, and lower bins in disk size offset ARd, 
respectively). Best-fit linear relations for each disk size 
bin are shown by colored dashed lines and best-fit pa- 
rameters are listed in Table |16[ The slopes of the re- 
lations are less than unity for both stellar mass esti- 
mates: (Mdyii/M«)opt decreases with stellar mass. In 
subsequent figures, we will show only results for M,,Beii. 

Figure 31 shows the relation between (Mdyn/A/,)opt 



(top panel) and (Mdyn/Mbar)opt (bottom panel) and 
stellar masses Mt,opt. Here, different symbols show bins 
in disk size offsets A_Rd. Figure[32]shows the correspond- 
ing relations with Rgg. Here, different symbols show bins 
in stellar mass Mt,Beii. In both figures, red triangles, yel- 
low circles, blue inverted triangles show the upper, mid- 
dle, and lower bins, respectively. Best-fit relations are 
shown by dashed lines of corresponding colors in Fig. [3l] 
and the best-fit parameters are listed in Table [16] 

We begin by discussing the trends in (Mdyn/M,)opt 
(top panels of Figs. 31 and 32 I. Then, we discuss the 
trends in (Afdyn/Afbar)opt (bottom panels of the same 
Figures). Our interpretation of these trends will be dis- 
cussed in §12. 4| 

We confirm that at a fixed stellar mass, 
larger (smaller) disks tend to have higher (lower) 
(Mdyn/M.)opt ( jZavala et al.|2003| [Pliagno et al.|2005[ ). 
In addition, we find that while most disks show a clear 
trend in (A/dyn/Af*)opt with stellar mass, the small- 
est disks fall off that main relation and tend to have 
low (Afdyn/M,)opt regardless of their stellar mass (blue 
inverted triangles and dotted line in the top panel of 
Fig. 31 and leftmost points in top panel of Fig. 32 1 ^ 

For stellar mass estimates Af,_Bcii, with a fixed 
Kroupa IMF, intermediate and large disks have 
(Mdyn/M«)opt decreasing from « 10 to 3 as Af.^Bcii in- 
creases from 10^ to IO^^Mq; while the smallest disks 



This result is more apparent here than in Fig. 6 of |Piza-| 
|gno et af] ( |2005[ l because of the addition of newly-observed 
galaxies with low stellar masses and small disks in our sam- 
ple. 
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Figure 30. Relation between dynamical mass Mdyn,opt ^'Ud 
stellar mass Af» opt within the optical radius Rgo for 189 
galaxies in the child disk sample. Top and bottom panels 
show different stellar mass estimates M.^Bell and Mt^MPAi 
respectively. Different symbols show galaxies in three bins 
in disk size offsets A_Rd (red triangles, yellow circles, blue 
inverted triangles, for the upper, middle, and lower bins, re- 
spectively), with corresponding best-fit relations shown by 
colored dashed lines. Solid lines shows Afdyn,opt 



have (Afdyn/Af*)opt ~ 2.7, regardless of stellar mass[^ 
At the high M, end, the ratios converge to the same 
low value, regardless of disk size. In other words, if we 
ignore the contribution of the gas mass, maximal disks 



Since M, mpa is systematically lower than Af, Belli ex- 
cept at the highest mass end (c.f. Fig. [6] in i ]5.4.3[ l, the in- 
ferred (Afdyn/Af*)opt would be correspondingly higher. 
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Figure 31. Relation between stellar masses within Rgo, 
A/,, opt I and dynamical-to-stellar and dynamical-to-baryonic 
mass ratios within Rgg (top and bottom panels, respec- 
tively). Different symbols correspond to bins in disk size off- 
set Ai?d (red triangles, yellow circles, blue inverted triangles, 
for the upper, middle, and lower bins, respectively). 



Figure 32. Relation between Rgo and dynamical-to-stellar 
and dynamical-to-baryonic mass ratios within Rgo (top and 
bottom panels, respectively). Different symbols correspond 
to bins in stellar mass Mt Bcll (red triangles, yellow circles, 
blue inverted triangles, for the upper, middle, and lower bins, 
respectively). 



(i.e., with Afdyn,opt = Af,,opt) require a non-uniform in- 
crease in the IMF normalization varying from ~ 0.4 dex 
(for the highest M« and lowest Tid galaxies) to ~ 1 dex 
(for the lowest M« galaxies with intermediate-to-high 
i?d). 

Now, looking at (A/dyn/A/bar)opt, we find the same 
dependence on disk size, at a fixed M«,opt, seen in 
the corresponding relation for (Mdyn/M*)opt (bottom 
and top panels of Fig. 311. The relation between 
(Mdyn/Mbar)opt and Rso is actually tighter than the 
corresponding relation with (Mdyn/M«)opt (bottom and 



top panel of Fig. 32 1, because the trend with stellar mass 



(at fixed Rso) seen in the latter is not found in the for- 
mer. 

We find that (Mdyn/Afbar)opt is roughly constant 
with stellar mass, at a fixed ARd, unlike (Afdyn/A/, )opt. 
This is a direct consequence of the sign and size of the 
trend of increasing gas mass fraction as stellar mass de- 
creases (see Fig. [7|. For the full child disk sample, the 
mean value of (Afdyn/Afbar)opt ~ 2.6, which corresponds 
to a 44 per cent baryon contribution to the mass in- 
side Rso- Taking these baryonic masses at face value 
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(i.e., assuming that the gas and stars have the same 
spatial extent), maximal disks (i.e., with Mdyn,opt = 
Aibar.opt) require a non-uniform increase in the IMF 
normalization varying from ~ 0.4 dex (for the highest 
M, galaxies) to ~ 0.9 dex (for the lowest M, galax- 
ies); the amounts are larger for the low-mass galaxies 
because they have higher gas mass fractions and the 
photometrically-derived gas masses are unaffected by 
changes in the IMF. 

The scatter in (Mdyn/Mbar)opt at fixed M*,opt is 
smallest at the highest M« end, as is the case for 
(Mdyn/A'/*)opt. The measured scatter in (Mdyn/Mbar)opt 
is 0.23 dex for the fuU sample, and 0.22, 0.15, and 0.21 
dex for the bottom, middle, and top bins in disk size 



offsets, respectively (Table 16 1 



The correlation of both (Mdyn/M,)opt and 
(Mdyn/Mbar)opt wlth disk size at fixed stellar mass sug- 
gests that a combination of these two properties can 
lead to a tighter relation with Afdyn/Af*. We confirm 
that the stellar surface density within Tiopt, ^*,opt = 
A^.,opt/(27r7?opt) forms a tighter relation with Mdyn/M, 
than either M, or Ropt alone ( [Zavala et al.|[2003l |Piza-| 
gno et al.|2005| . 



The top panel of Figure |33] shows the relation be- 
tween E«,opt and (Mdyn/Af*)opt for the child disk sam- 
ple. The best-fit relation is shown by the dashed line 
and best-fit parameters are listed in Table [161 The bot- 
tom panel of Figure |33] shows the relation between 
the baryonic surface density within Rso, Sbar.opt = 



S.,opt X (Mbar/M.) and (Mdyn/Mbar)opt. As in Fig. 31 



different symbols correspond to bins in disk size offsets 
ARd (red triangles, yellow circles, blue inverted trian- 
gles show upper, middle, and lower bins, respectively). 

Comparing the upper panels of Figs. |31| and|33[ we 
find that the smallest disks, which fall below the relation 
involving M, move to the right (toward higher E*,opt) 
by an amount that places it roughly on the Et,opt- 
(Mdyn/M«)opt relation defined by the larger disks. 



•,opt (Mdyn/Mbar)opt is COm- 
■,opt~ 



The scatter in the Eba 
parable to the scatter in the S»,opt-(A/dyn/Af,)opt 
relation — equal to 0.18 and 0.15 dex, respectively. The 
trends in (Mdyn/Mbar)opt with M,,opt and Rso (lower 



panels of Figs. 31 and 32 \ indicate that the tightness of 



its relation with Ebar,opt is almost entirely driven by its 
dependence on disk size, rather than stellar mass. This 
observation has important consequences for our inter- 
pretation of these trends in §12. 4| 

Finally, to look for the tightest possible relation 
for (Mdyn/Af*)opt, we explore other combinations of the 
form Af,, opt /-Ropt) with the exponent k treated as an ad- 
ditional free parameter. Maximum likelihood fits yield 
K — 1.7 ± 0.3, consistent with the definition of the stel- 
lar surface density (with n — 2). Therefore, we conclude 
that the stellar surface density is, indeed, a near-optimal 
photometric tracer of (Mdyn/Mt)opt. 
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Figure 33. Top panel: Relation between dynamical-to- 
stellar mass ratio (Afdyn/Mt)opt and stellar surface density 
within the optical radius ijgo, S.^opt- Bottom panel; Relation 
between dynamical-to-baryonic mass ratio (A'/dyn/A'/bar)opt 
and baryonic surface density within the optical radius -RgOi 
Sbar,opt- In both panels, red triangles, yellow circles, blue 
inverted triangles show upper, middle, and lower bins in disk 
size offset A_Rd , respectively. The dashed line shows the best- 
fit linear relation for the full sample and dot-dashed lines are 
displaced from this relation by ±5- to show the amount of 
intrinsic scatter. 



12 SUMMARY AND DISCUSSION 
12.1 Methodology 

Our sample selection and methodology are dictated by 
our primary goal of obtaining photometric estimates of 
disk rotation velocity for a large sample of galaxies with 
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Table 16. Fits to the M^y^ — Aft^opt relation and M^yn/Af, relations with various physical parameters. 
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0.174 










middle 1 /3 in AR^ 


10.07 


10.676 (0.010) 


0.783 (0.030) 


0.062 (0.020) 


0.114 










upper 1/3 in A/?d 


10.07 


lU.ozz (U.OzO) 


O.TOo (0.044) 


n lot" / n r\ooN 

0.125 (0.028) 


0.151 


log 




log 


M«,MPA,opt 


all galaxies 


O OQ 


1 n (r\ m 
iU.DDl t^U.Uio^ 


u.oou yyj.vzy ) 


U. toy l^U.UiD j 


n ono 
v.zyjz 










lower 1/3 in A/2d 


y .yo 


in /ini noQ\ 
iU.4y 1 [d.VZo ) 


U. / 40 (^U.U04 ) 


n 1 on fn noQ\ 
U. Izy (^U.Uzo j 


n 1 on 
U.loU 










miaule i/ o m ZAitd 


y .yo 


iU.OOy t^U.Ulo ) 


U.D ( / (^U.UoU ) 


U.UUU (^U.Uzy J 


U.llD 










upper 1/3 in AJ2d 


y .yo 




U.OOO l^U.UOO ) 


n 1 9*? f'n ^^^^ 

U. IZO l^U.UOU ) 


U. 141 


log 


Mdyn/Af* 


log 


A^,, Bell, opt 


all galaxies 


10.07 


0.600 (0.014) 


-0.225 (0.037) 


0.168 (0.015) 


0.202 










lower 1/3 in A_Rd 


10.07 


0.428 (0.023) 


-0.080(0.061) 


0.136(0.019) 


0.174 










middle 1/3 in ARd 


10.07 


0.608 (0.011) 


-0.215(0.030) 


0.057(0.028) 


0.114 










upper 1/3 in AiJd 


10.07 


0.753 (0.018) 


-0.309 (0.040) 


0.113(0.028) 


0.150 


log 


Afdyn/Af* 


S, 


opt 


all galaxies 


7.33 


0.603 (0.009) 


-0.419(0.027) 


0.097(0.012) 


0.145 


log 


Mdyn/Afbar 


log 


Boll, opt 


all galaxies 


10.07 


0.358 (0.017) 


-0.004(0.045) 


0.221 (0.023) 


0.232 










lower 1/3 in A_Rd 


10.07 


0.199 (0.026) 


0.170(0.084) 


0.203(0.035) 


0.215 










middle 1/3 in A_Rd 


10.07 


0.362(0.017) 


-0.011 (0.041) 


0.127(0.018) 


0.154 










upper 1/3 in Ai?d 


10.07 


0.499 (0.026) 


-0.098 (0.060) 


0.190(0.029) 


0.206 


log 


Mdyn/Afbar 


log 


^bar.opt 


all galaxies 


7.33 


0.470 (0.014) 


-0.432(0.045) 


0.140(0.013) 


0.177 



Notes. — Col. (1): y is the dependent variable in the fit. Col. (2): x is the independent variable in the fit. Col. (3): Subsample 
or bin used for the fit. Col. (4): Xp is the pivot value for x. Cols. (5-7): Best-fit parameters and their l-cr uncertainties: a is the 
zero-point, b is the slope, and a is the intrinsic Gaussian scatter. Col. (8): a"meaa is the measured scatter in the mean relation, 
defined to be the RMS of the residuals (Aj/); = yi — [a + b X {x — Xp)]. Masses are in units of Mq, and E^^opt and Sbar,opt ^^e 
in units of Mq kpc"'^. 



imaging data. To achieve this goal, we construct a TFR 
sample of 189 galaxies that is a fair subsample of a well- 
defined parent disk sample of ~170 000 galaxies from 
the SDSS (!|3|. 

To derive robust TFRs, we use rotation velocity 
amplitudes Vso defined at Rso, the radius containing 80 
per cent of the i-band galaxy light and approximately 
at the peak of the total (disk-|-DM halo) rotation curve 
for typical disk galaxies and their haloes (c.f. S 6.5.3 and 
S 9.1 1. We calculate and account for systematic errors in 



Vgo, in addition to statistical errors, to downweight the 
contribution of those few galaxies with peculiar rotation 
curves to the TFR fits. We also perform weighted fits to 
the TFR to account for the small difference between the 
stellar mass functions of the TFR and parent disk sam- 
ples; however, in practice, the two are similar enough 
that the weights do not significantly affect the fits. 



12.2 Tests of systematics 

We find that the TFR is robust to many systematic ef- 
fects, including slit misalignments (S 6.5.1 1, differences 
in observing instruments and conditions 6.5.2 1, dif- 
ferences in analysis pipelines (i.e., the rotation curve 
extraction and fitting methods, and bulge-disk decom- 
position fits, which differ between this work and P07) 



(i 6.5.3 and ^9.51, and differences in sample selection 
(i.e., between this work and P07) (|9.2| and j ]9.5| |. We 
also verify that our internal extinction corrections suc- 
cessfully remove the weak correlation between disk axis 



ratios gd and velocity residuals from the (uncorrected) 
Mf^ ITFR (j ]l0.l| . This indicates that even if they are 
uncertain on a galaxy-per-galaxy basis, they are ade- 
quate for the TFR sample on average. 

While there are many other potential sources of sys- 
tematics that we have not explicitly tested for, we note 
that the observed small scatter in the TFR — 0.056 dex 
of total scatter in log Vgo, corresponding to 0.20 dex in 
log A'/«, Boll or 0.5 mag in Mi — itself serves as an upper 
bound for their contribution. We find that the expected 
additional sources of scatter (e.g., intrinsic disk elliptici- 
ties) can explain essentially all of the observed scatter in 
the TFR, regardless of whether our measurement errors 



were overestimated (^7.31. The disk-dominated sample 
yields an even smaller total TFR scatter of 0.05 dex in 
log Vgo, corresponding to 0.18 dex in logAf.^Boii (or 0.44 
in magnitude "units"). For comparison, previous TFR 
studies based on pruned samples of spirals (for distance 
indicator work) found a range of observed total scatter 
~ 0.3-0.4 mag ( Strauss fc Willick||1995 and references 
therein). 



12.3 Optimal estimator of disk rotation 
velocity 

One of the primary aims of this work is to identify and 
calibrate an optimal estimator of disk rotation velocity, 
i.e., one with minimal scatter in its ITFR. We confirm 
that the use of redder bands C Strauss fc Wimck||1995|), 
or the addition of colour information ( Kannappan et al. 
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2002 1, yield tighter ITFRs. It is interesting to note that 



A^syn(A, (? — r) for different bands yielded not only IT- 
FRs with identical scatter, but almost identical slopes 
and zero-points as well. We note that Msyn can be inter- 
preted as an extrapolation to some redder (i.e., infrared) 
band, which traces stellar mass better than any single 
optical band. We also find that the best-fit coefficients 
a for the colour terms in Msyn{^,g — r) are consistent 
with the corresponding coefficients a^'^" from the Mt/L 
fitting formulae of Bell et al. (2003) (see Table[9|. Then, 
in §10. 2| we explicitly show that almost all of the corre- 
lation between g — r colour and velocity residuals from 
the Mi ITFR can be attributed to, and hence removed 
by, the colour dependence of the stellar mass-to-light 
ratio, M./Li oc (g - r)° ®'^ (also from BeU et al. 2003). 

Ultimately, we choose Af*,BGii to be our optimal es- 
timator since it has a natural physical interpretation 
(unlike Af^yn) and is straightforward to calculate (un- 



like M*,mpa). We find (c.f. Table 10 and Fig. 201 



log V8o(Af,,Bcii) = (2.142 ± 0.004) + (0.278 ± 0.010) 
X (log A/.,Beii - 10.102), (39) 

with an intrinsic Gaussian scatter a = 0.036 ± 0.005 
in log Vgo (using a Kroupa IMF to define Af*^Beii). The 
associated uncertainty in this estimate is given by the 
measured scatter in the ITFR, cTmeas = 0.056 dex. 

Pruning the galaxy sample to include only disk- 
dominated galaxies is another way to obtain a tighter 
ITFR. We show in §9.4| that removing galaxies with 
D/T ^ 0.9 yields a tighter Af.,Beii ITFR. We note that 
pruning a large imaging dataset with many fainter and 
smaller galaxies than those in our sample is nontrivial, 
but worth investigating. 



12.4 Relative contributions of stars and DM 

The question of the relative contribution of stars and 
DM in the optical regions of disk galaxies remains unre- 
solved. Observed rotation curves only weakly constrain 
the relative contributions of the different components: 
stellar and gas disks, bulge, and DM halo. To disentan- 
gle these components, we assume a fixed stellar IMF, 
then infer stellar and gas masses based on the galaxy 
luminosities and colours. We use our observations to 
confirm several known trends and point out interesting 
new ones that will help constrain the possible scenar- 
ios. In future work, we will aim to construct disk galaxy 
models that reproduce these observations. 

In §10. 3| we study correlations between velocity 
residuals from the Af« ITFR and disk size offsets A_Rd 
for the child disk sample divided into three bins in Af, 
(Figure 291. We find that the lowest stellar mass galax- 



ies show a positive correlation between TFR residuals 
and disk size offsets, suggesting that they have a domi- 
nant NFW DM halo. Moreover, the sign of the correla- 
tion changes from positive to negative toward increasing 
stellar mass, suggesting a progression toward an increas- 
ing contribution of the stellar disk. (Recall that a pure 



NFW DM halo yields a slope of +0.5, while a pure ex- 
ponential disk yields a slope of —0.5.) 

In §11[ we use stellar mass estimates Aft^Beii and 
baryonic mass estimates Afbar (defined in [ 5.4 and [ 5.5 



respectively) with a fixed Kroupa IMF, to calculate 
dynamical-to-stellar and dynamical-to-baryonic mass 
ratios within Rso, (Mdyn/Af*)opt and (Afdyn/Afbar)opt, 
respectively. We find that (Afdyn/Af*)opt decreases from 
~ 10 to 3 as A'/,_Beii increases from 10^ to W^^Mq (top 



panel of Fig. 311. This trend of an increasing (decreas- 
ing) contribution of stars (DM) to the potential on op- 
tical scales as stellar mass increases is consistent with 
that indicated by the size residual correlations. 

On the other hand, because low M, galaxies 
have higher gas mass fractions than high Af, galaxies, 
(A/dyn/Afbar)opt is roughly constant with stellar mass, 
with a mean of 2.6 and standard deviation of 0.23 dex 



(bottom panel of Fig. 311. Blanton et al. (20081 found 



a similar result from a sample of isolated dwarf galax- 
ies with luminosities M,. — 51og/i = —14 to —21.5 mag. 
At face value, it seems that this finding is inconsistent 
with the observed dependence of the sign of the size 



residual correlations on stellar mass (f 10.3 1. However 



recall that we have not taken into account the differ- 
ence in the spatial distributions of the gas and the stars. 
In fact, (Afdyn /Afbar) opt and (A/'dyn/Af*)opt bracket the 
more realistic scenario, in which only some fraction of 
the gas mass is within the optical radius. We defer de- 
tailed analysis and discussion to a future modelling pa- 
per. 

Panels of Figure |3l] clearly show that the scatter in 
both (A/dyn/A/*)opt and (A/dyn/Afbar)opt at fixed stellar 
mass is driven by differences in disk sizes: larger disks 
tend to have higher values in both ratios. As a conse- 
quence, both ratios form tighter relations with stellar (or 
baryonic) surface density than with stellar mass or disk 
size alone ( Zavala et al. 2003 ; Pizagno et al. 2005 ) . In the 



disk galaxy models of jFirmani fc Avila-Reese^(2000J , the 
observed tight trend of (Afdyn/Af,)opt with Et,opt arises 
from a correlation of star formation efficiency with sur- 
face densityp^ Here , we propose an alternative explana- 
tion: the observed variation in (Afdyn/Af,)opt with disk 
size is actually tracing the slope of the (average) to- 
tal mass profile of disk galaxies at these radii. In our 
own disk galaxy models, we will show that the observed 
trends are reproduced even if the disk mass fraction is 
independent of disk size or stellar surface density (R. 
Reyes, in prep.). 

The reason this happens, qualitatively, is that if we 
assume a fixed DM halo for a given disk (ignoring real- 
istic variation due to, e.g., scatter in halo concentration 
at fixed mass) , then the larger disks will necessarily have 
more DM within Rgo . This assumption is not necessarily 
true- it is possible, for example, that a disk that is larger 

^® This has, in fact, motivated disk galaxy models in which 
the disk mass fraction (the ratio of the stellar mass to the 
total halo mass) is dependent on S,,opt (e.g., [Gnedin et al.| 
20071. 
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or smaller than the typical size for its stellar mass would 
have a different effect on the DM halo density profile. 
Our results suggest that such a correlation between DM 
halo profile modifications and disk size at fixed M« are 
not of the right sign, or sufficient magnitude, to cause a 
deviation from our naive expectation that the disk size 
simply modulates how much of the dark matter the gas 
in the disk "sees." 



12.5 Applications to future work 

The calibrated TFRs presented here can be used for a 
variety of purposes: for estimating disk rotation veloci- 
ties for large photometric samples, connecting galaxies 
with their DM haloes, and studying models of galaxy 
formation and evolution. As mentioned above, we will 
use these observations to construct and constrain mod- 
els of disk galaxy formation. In another paper, we will 
constrain Kot/V^ir as a function of stellar mass, using 
our calibrated ITFR to determine Vi-ot , and weak lensing 
measurements to determine Vvir- As discussed in [jT] this 
measurement constrains the shape of the total mass pro- 
file, and can potentially indicate whether baryons have 
modified the gravitational potential well, in which sense, 
and by how much. 
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APPENDIX 

In this Appendix, we derive the relation between 
dynamical-to-stellar mass ratio within a radius R 
(A/dyn/M*)(_R) and the observed rotation velocity at 
that radius Kot (Eq. 361. First, we consider an ideal- 
ized galaxy model with two components: an exponen- 
tial stellar disk and a NEW DM halo (we neglect the 
bulge for now, for simplicity). The disk has an infinites- 
imal thickness, scale length i?d, central surface mass 
density Eq, and a surface mass density distribution, 
Ed(-R) = Eo exp(— _R/i?d). The mass interior to R is 
given by 



Md(i?) 



2tvT,oRI [1 ~ exp(-i?/i?d) (1 + R/Rd)] 
2nY;oRlT{R/Rd). (40) 



and the rotation curve Vc,d{R) is given by (Freeman 
T970| ) 

4^GEoi?dy'[/o(y)i^o(y) - /i(y)^i(y)] 
47rGEoi?de(y), (41) 



where y = R/{2Rd), G = 4.3012 x 10"*^ is the gravi- 
tational constant in units of kpc (km/s)'^ Mq^, and li 
and Ki are the modified Bessel functions of the first 
and second kind. The disk rotation curve will peak at 
-Rdisk.max = 2.15i?d and the peak velocity is about 15% 
higher than that of a spherical mass distribution with 
the same interior mass. 

The NEW halo has an interior mass Mnfw(-R) 
and a rotation curve Vc,nfw{R)- If one (incorrectly) 
assumed that the observed (i.e., total) rotation veloc- 
ity Kot arises from a spherical mass distribution, one 
would infer a dynamical mass Afsph interior to R given 
by 



Mspb{R) 



.2 +v;2^) 



RV^R) _ R{Vl 



G G (^2) 

Mt,t,w{R) + K.{R/Rd)Md{R), (43) 



where we have introduced a geometric correction factor 



ICiR/Rd) 



Vc.d(R) 
K'd,sph(-R) 



47rGEofldg(y) 
GMd{R)/R 

47rGEofldg(2/) 
2TTGT,oRlT{R/Rd)/R 
2B{y) R 



T{R/Rd) Rd ' 



(44) 
(45) 
(46) 



Here, 14,d,sph is the rotation curve for a spherical mass 
distribution with the same interior mass as the exponen- 
tial disk, and T and B are functions defined in Eqs. |40| 
and |41| respectively. We find that K. peaks at around 
~ 3i?d, with a value of 1.3, and then falls off slowly 
and asymptotes to 1. For R/Rd = 2.2, {F,B,1C) = 
(0.645,0.1935,1.319); for R/Rd = 3, the values are 
(0.8009,0.1795,1.345). 
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The dynamical-to-steUar mass ratio is therefore re- 
lated to its "spherical" counterpart by 

Mdyn(i?) _ , , MNFw(i?) 



Ma{R) 



M,pi,{R) 
Mi{R) 



1 - K. 



(47) 
(48) 



In other words, if one assumes a spherical mass distri- 
bution, one would falsely attribute a larger dynamical 
mass to the stellar disk, and therefore overestimate the 
total-to-stellar mass ratio by a constant additive term 
K.{R/Rd) - 1 « 0.3 for R/Ra = 2-3. 

We now generalize the above results for a model 
galaxy with an additional bulge component. Since the 
bulge mass distribution is approximately spherical, it 
can be lumped together with the halo mass in Eq. |43[ 
Defining Afstar = Md + Mbuigo, we have 



Mdyn(-R) 
Mstar(-R) 



[MNFw(fi) + Mbuigc(-R)] + MdjR) 

Metar(-R) 

Msph(fl) 
Mstar(i?) 

MjjR) 
Mst^.iR) 



K. 



R^ 
Rd 



(49) 



In practice, we adopt the radius R — Rso and 
Md(-R8o)/Matar(-R8o) = {D/T)m, the disk-to-total light 
ratio interior to R^o (Eq. |37[ ). The correction is sup- 
pressed in galaxies with significant bulges, because the 
contribution of the disk is smaller than what was as- 
sumed when all the light was attributed to the disk. 
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